August 30, 2010 0:10 WSPC/INSTRUCTION FILE electron-positron 



> 



International Journal of Modern Physics D 
© World Scientific Publishing Company 



o 

(N 

^ ■ Excesses in the Cosmic Ray Spectrum and Possible Interpretations 

< 

Yi-Zhong Fan 

Department of Physics and Astronomy, University of Nevada, Las Vegas, NV 89119, USA; 
Purple Mountain Observatory, Chinese Academy of Sciences, 210008, Nanjing, China 

X 

Oh; 

Q . Jin Chang 

^ Purple Mountain Observatory, Chinese Academy of Sciences, 210008, Nanjing, China 

C/3 ' 
. 

Received Day Month Year 
Revised Day Month Year 
Communicated by Managing Editor 



Bing Zhang 

Department of Physics and Astronomy, University of Nevada, Las Vegas, NV 89119, USA 



The data collected by ATIC, PPB-BETS, FERMI-LAT and HESS all indicate that 
there is an electron/positron excess in the cosmic ray energy spectrum above ~ 100 
GeV, although different instrumental teams do not agree on the detailed spectral shape. 
PAMELA also reported a clear excess feature of the positron fraction above several GeV, 
00 ' but no excess in anti-protons. Here we review the observational status and theoretical 

C " 3 ' models of this interesting observational feature. We pay special attention to various 

C " 3 , physical interpretations proposed in the literature, including modified supernova rem- 

nant models for the e^ background, new astrophysical sources, and new physics (the 
dark matter models). We suggest that although most models can make a case to inter- 
pret the data, with the current observational constraints the dark matter interpretations, 
k>( ^ especially those invoking annihilation, require much more exotic assumptions than some 

. astrophysical interpretations. Future observations may present some "smoking-gun" ob- 

' servational tests to differentiate among different models and to identify the correct in- 

terpretation to the phenomenon. 

Keywords: excess. Supernova remnants, Pulsar, Dark matter. Particle 

1. Introduction 

Cosmic rays (CRs) are energetic particles originating from outer space that impinge 
on Earth's atmosphere. Almost 90% of all the incoming cosmic ray particles are 
protons, about 10% are helium nuclei (alpha particles), and slightly under 1% are 
heavier elements and electrons^. The history of the CR studies may be rooted in 
1901 when two groups came to the same conclusion that pure air in a closed vessel 
possessed some electrical conductivity^^'^, while no visible sources of air ionization 
were observed. At that time it was already known that X-rays and radioactivity were 
factors contributing to the enhanced electrical conductivity of gases. Therefore, the 
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observed effect of "dark current" caused by residual ionization of the air, was then 
regarded as being associated with radioactive contamination both in the air and in 
the environment^. Wilson speculated that the residual ionization was due to (x^rtain 
highly penetrating radiation coming from outside the Earth's atmosphere*. Such a 
speculation was convincingly confirmed by Hess^ and Kolhorster^ about 10 years 
later. August 7, 1912, the date of Hess' most successful balloon flight among the 
total of ten, is generally regarded as the date of the discovery of cosmic rays. The 
discovery of cosmic rays opens a new window not only for high energy astrophysics 
but also for high energy physics. This is because prior to the constructions of modern 
particle accelerators, high energy particles could only be observed in cosmic rays. 
Between 1932 and 1953, many new particles, including positrons, /i^-leptons, tt^- 
mesons, K^— and /sT"— mesons, A°— , S~— and S+— hyperons, were all discovered 
in use of cosmic rays (Ref.l and the references therein). 

A thorough review of CRs can be found in many excellent articles and 
books'''^'^'^°, and is beyond the scope of this review. Here we will focus on a small 
fraction of these high energy particles, the electron/positron (e^) CRs. Although 
CRs had been discovered in 1912, CR electrons reaching the earth at energies above 
a few hundred MeV had not been convincingly identified until 1961^^'^^. The first 
positron-to-electron ratio in CRs up to an energy ~ 1 GeV was later reported in 
1964^^. Most CR electrons are likely from the supernova remnants while the CR 
positrons are mainly produced through hadronic processes when CR protons collide 
with intergalactic hydrogen. In siich a scenario there should be no prominent feature 
in the TeV energies in the electron/positron total spectrum, and the positron-to- 
electron ratio TZ = + $e-) should drop with energy monotonously^*. A 
rise of the positron-to-clectron ratio above ^ 10 GcV was observed in 1974'^^''. Such 
a tendency was confirmed by several balloon flights, including the 1976 flight by 
the New Mexico State University group launched from Palestine, Texa in 1976^°, 
and 1984 flight by the University of Chicago group launched from Hawaii^^, and in 
particular, by the recent PAMELA mission^^^^^ (see section 2 for extensive discus- 
sion). In 1987 it was already clear that the rising behavior of TZ above ~ 10 GeV 
is not predicted by calculations of CR propagation^^'^^ and may imply an excess 
of CR positrons. Recently, ATIC^*, PPB-BETS^^ and Fermi-LAT^s reported elec- 
tron/positron total spectrum up to TeV and found a hardening/bump in the energy 
range from 100 GeV to 1 TeV. In the conventional approach, the injection spectrum 
of the electrons (positrons) is taken as a single power-law. Since diffusion and elec- 
tron/positron cooling are more efficient in higher energies, one would expect that the 
spectrum should soften with energy. Therefore the unexpected spectrum hardening 
observed by these detectors strongly suggest another excess (besides the positron 



*In 1965, a group reported TZ ~ 0.8 at 10 — 30 GeV, implying an excess of positron CRs^^. It was 
however only based on a total of 13 events and therefore highly uncertain ^'^ . If one plots the single 
data of Ref.18 (at ~ 10 GeV) published in 1969 together with some other data at lower energies 
reported before, one would also get a positron excess^®. 
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fraction excess). Interestingly this total excess is consistent with extrapolating 
the positron fraction excess to high energies. These interesting features draw a lot 
of attention, and various physical origins have been explored. In the literature, new 

astrophysical sources or even new physics (dark matter) are widely suggested. Some 
more conservative authors argue that the excesses are simply due to inadequate ac- 
counts on the electron/positron CR background in previous modeling. The situation 
is unclear. Several reviews appeared recently but are mainly focused on dark matter 
models ^7,28 ^^lis review we plan to present a more complete overview of this 
interesting observational feature and various physical interpretations. 

The structure of this review is as follows. We first discuss the observational 
aspects of the excesses in Section 2. We then review the cosmic ray propagation 
models and the general cosmic sources of in Section 3. In the next three sections 
(4-6), we discuss three proposed physical origins of the observed excess, including 
the modified supernova remnant (SNR) models of background, new astrophysical 
sources, and new physics. We summarize the strengths/weaknesses of various models 
in Section 7 with a commentary on the prospects. 

2. Observations 

Different from the CR protons/nuclei, the TeV electrons/positrons suffer significant 

energy loss and cannot reach us if the source is not nearby (i.e., < 1 kpc). Accu- 
rate measurements of the high-energy electron/positron CRs then provide a unique 
opportunity to probe the origin and propagation of CRs in the local interstellar 
medium and to constrain models of the diffuse gamma-ray emission. Observations 
of electrons have been notoriously difficult because of their low intensity requiring, 
in particular at high energies, rather large detectors, and because of the need of 
effective discrimination against proton-induced background. So far the high-energy 
electron spectrum was measured by some balloon-borne experiments (in addition to 
those already mentioned in the introduction, some others can be found in Refs. 29, 
30, 31, 32, 33, 34, 35, 36, 37, 38), some space missions (including IMP-l^s, IMP-III^o, 
IMP-IV^i, IMP-7*2, AMS-01^^ PAMELA^i and Fermi^^), and the ground-based 
Cherenkov telescope H.E.S.S.^*. 

2.1. Experiments 

2.1.1. Balloon-borne experiments 

We first briefly review the development of technology used in the experiments. In 
two of the very early Balloon-borne experiments^^'^^, the cloud chamber technol- 
ogy was employed, which was soon abandoned with the introduction of other more 
sophisticated detectors, both visible and electronic (see Ref.l7 for a review on the 
development of technology in the early experiments). In the two flights at Fort 
Churchill in 1963^^, the spark chamber technique was adopted. With the help of 
one permanent magnet, electrons and positrons were distinguished. Such an instru- 
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ment was later improved^ ^, through replacing the liquid Cerenkov detector by a 
gas Cerenkov detector to enable the determination of the momentum and charge of 
electrons up to 5 — 10 GeV. A group from Bombay^^'*^ adopted the nuclear emulsion 
stack technique, which measures CR electrons to an energy above 50 GeV. In order 
to identify electron showers efficiently, multiple pulse-height analyses of many coun- 
ters and time of flight techniques had been used in 3 high-latitude balloon flights 
from Palestine, Texas in 1970 by the Chicago group^^. Later they added a transi- 
tion radiation detector, consisting of a sandwich of six radiators and six multiwire 
proportional chambers, to improve electron and positron discrimination^^'^^. In a 
balloon-borne magnetic spectrometer experiment^^'^^, a super-conducting magnetic 
spectrometer was adopted, and a new bremsstrahlung-identification technique was 
developed to measure the separate e~ and e+ spectrum in the primary CRs from 
4 to 50 GeV. The emulsion chamber technique was used to detect CR electrons 
to an energy above 1 TeV^"'^^. In the experiment of Balloon-borne Electron Tele- 
scope with Scintillating Fibers (BETS), the absolute energy spectrum of electrons 
was measured with a highly granulated fiber calorimeter^*. Finally the Advanced 
Thin Ionization Calorimeter (ATIC) experiments used an ionization calorimeter to 
measure the energy deposited by cascades formed by particles interacting in a thick 
carbon target 

Since 1950, dozens of balloon-born experiments dedicated to electron/positron 
CR observations have been carried out. Below we only mention some recent ones 
that arc most relevant to the main topic of this review, the electron/positron excess. 

Two early charge determining experiments. The Chicago group firstly 
used a permanent magnet to measure the positron-to-electron ratio. In the energy 
range 0.2 GeV < 5 < 10 GeV, the ratio TZ was found to decrease with energy"^^. 
Adopting a super-conducting magnetic spectrometer, TZ was measured from 5 GeV 
to 50 GeV^^. Although scatter is large, an tendency of increasing TZior £ > 10 GeV 
is present. Such a behavior was confirmed by some subsequent groups (e.g. Ref. 19, 
20), but not by some others. 

The emulation chamber experiments by a Japanese group started in 
1968. The e+ -|- e~ spectrum derived by combining data from the cham- 
bers exposed from 1968 to 1976 was well represented by J{E) = 1.6 x 
10-''(£:/100 GeV)-3-3±"-2 sr"! s'^ GcV^i in the energy range of 30 - 10^ 
GeV^°, which is steeper than the spectrum J cx £:-2.6i±o.i obtained in the early nu- 
clear emulation stack experiment by the Bombay group^^. This confirms some early 
results, for example, obtained with an ionization calorimeter'''''^'''^. The observations 
in 1996 and 1998 extended the e"*" -|- e~ spectrum to ~ 2 TeV^^ and might suggest 
a steepening above 1 TeV. 

High-Energy Antimatter Telescope (HEAT) is a NASA-supported pro- 
gram of high-altitude balloon-borne experiments to study antimatter in the pri- 
mary cosmic radiation. The HEAT detector^® consists of a magnetic spectrometer 
combined with a transition radiation detector, an electromagnetic calorimeter, and 
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time-of-flight scintillators. It can measure the cosmic-ray positron fraction from 1 
to 50 GeV. A new version of the HEAT instrument, HEAT-pbar was designed to 
observe the high-energy cosmic-ray antiproton flux but it is also suited for the ob- 
servation of electrons and positrons at energies below 15 GcV'^"''. The combined data 
from the three HEAT flights indicate a small positron flux of nonstandard origin 
above 5 GeV^^. The evidence for the tendency of increasing 7?. for f > 10 GeV is 
however rather weak, consistent with what was found in the Cosmic AntiParticle 
Ring Imaging Cherenkov Experiment in 1998^^. 

BETS and PPB-BETS. BETS represents a detector, the Balloon-borne Elec- 
tron Telescope with Scintillating fibers, which preserves the superior qualities of 
both electronic detectors and emulsion chambers With such a balloon-borne 
pay load, cosmic-ray electrons were observed in the energy range from 12 to ~100 
GeV. The energy spectrum is described by a power-law index of —3.00 ± 0.09, and 
the absolute differential intensity at 10 GeV is 0.199±0.015 m'^ s'^ sr"! GeV"^^^. 
PPB-BETS is an advanced version of BETS measuring the electron spectrum up 
to TeV and was flown by Polar Patrol Balloon (PPB) in Antarctica. The results 
from the PPB-BETS calorimeter indicates a possible hump in the energy range 100 
GeV to 1 TeV in the total spectrum^^, consistent with that reported by ATIC (see 
below). 

The Advanced Thin Ionization Calorimeter (ATIC) is a balloon-borne 

instrument flying in the stratosphere over Antarctica to measure the energy and 
composition of cosmic rays. ATIC was laimchcd from the McMurdo Station for the 
first time in December 2000 and has since completed three successful flights out of 
four. The detector uses the principle of ionization calorimetry: several layers of the 
scintillator bismuth gcrmanate emit light as they arc struck by particles, allowing to 
calculate the particles' energy. A silicon matrix is used to determine the particles' 
electrical charge^^. The main discovery is a spike-like e+ -|- e~ spectrum feature in 
the energy range of 300 - 800 GeV^^. 

2.1.2. Space station and satellite experiments 

The series of Interplanetary Monitoring Platform (IMP) experiments measure the 
low energy electron CRs {£ < 1 GeV), which are not the focus of this review. In 
the following we only focus on several relevant missions. 

AMS-01 (Alpha Magnetic Spectrometer - 01) is the predecessor to AMS-02, 
a detector to be operated on the International Space Station (ISS) for at least 
3 years (http://ams.cern.ch/). The AMS-01 experiment was flown on the Space 
Shuttle Discovery for ten days in June, 1998. It was designed to search for various 
types of unusual matter by measuring cosmic rays. The detector consisted of a per- 
manent Nd-Fe-B magnet, six silicon tracker planes, an anticoincidence scintillator 
counter system, the time-of-flight (TOF) system consisting of four layers of scintil- 
lator counters and a threshold aerogel Cerenkov detector. With these instruments, 
the high-energy proton, electron, positron, helium, antiproton and deuterium spec- 
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tra were accurately measured. The obtained TZ in the energy range of 1 — 30 GeV 
is largely consistent with that obtained by HEAT 

PAMELA, the Payload for Antimatter Matter Exploration and Light-nuclei 
Astrophysics, is a powerful particle idcntifior using a permanent magnet spectrome- 
ter with a variety of specialized detectors (http:/ /pamela. roma2.infn.it/index.php). 
It is measuring with unprecedented precision and sensitivity the abundance and en- 
ergy spectra of CR electrons, positrons, antiprotons and light nuclei over a very large 
range of energy from 50 MeV to hundreds GeV, depending on the species^^'^''. The 
most exciting discovery made so far is the unambiguous detection of the increasing 
behavior of TZ in the energy range 10 — 100 GeV 21.22 

Fermi 7— ray Space Telescope is an international and multi- agency space 
mission that studies the cosmos in the energy range 10 keV - 300 GeV 
(http://fermi.gsfc.nasa.gov/). The Large Area Telescope (LAT) is the main instru- 
ment on-board. Conceived as a multipurpose observatory to survey the variable 
gamma-ray sky between 20 MeV and 300 GeV including the largely unexplored 
energy window above 10 GeV, it is designed as a low aspect ratio, large area pair 
conversion telescope to maximize its field of view and effective area. The LAT an- 
gular, energy and timing resolutions rely on modern solid state detectors and elec- 
tronics. Since electromagnetic cascades are germane to both electron and photon 
interactions in matter, the LAT is also by its nature a detector for electrons and 
positrons^^. The measurement of the CR electron spectrum from 20 GeV to 1 TcV 
based on the first half year data also displays a spectrum hump in the 100 — 1000 
GeV range, although the signature is not as prominent as that found by ATIC and 
PPB-BETS. 

2.1.3. Ground-based telescopes 

H.E.S.S, the High Energy Stereoscopic System, is an array of four imaging at- 
mospheric Chcrcnkov telescopes located in the Khomas Highland of Namibia 
(http://www.mpi-hd.mpg.de/hfm/HESS/). It is designed to obtain the Cherenkov 
images of cosmic-ray hadrons and electrons, as well as gamma-ray photons (sorted 
according to the relative flux). To measure the spectrum of electrons, the other two 
components (hereafter background events) have to be eliminated. First, data taken 
from the sky regions known to contain gamma-ray sources are rejected. However, the 
most critical aspect of the electron analysis is the efficient rejection of the hadronic 
background. A Random Forest approach was used to convert the image information 
into a single parameter describing the degree to which a shower is electron-like, that 
spans from (for hardons) to 1 (for electrons). By Monte-Carlo simulations, the 
distribution of this parameter for simulated hadrons and electrons is determined 
(See Figure 1 in Aharonian et al.'''*). The number of measured electron showers in 
each energy band can then be deduced. In this way, the spectrum is measured from 
0.35 - 4 TeV. While the overall electron flux measured by H.E.S.S. is consistent 
with the ATIC data within statistical and systematic errors, the H.E.S.S. data do 
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not confirm a pronounced peak in the electron spectrum suggested by ATIC and 
PPB-BETS^^ 

2.2. Observation results 

Prior to 2008, the high energy e+ + e~ data are consistent with a featureless power- 
law spectrum within errors, which is in agreement with the theoretical predic- 
tions from both analytical and numerical models^"*. Recently, the ATIC balloon 
experiment^* discovered a prominent spectral feature at around 600 GeV in the 
6+ -I- e~ spectrum. Instead of a decaying feature expc;cted by electron/positron 
cooling, the team discovered a clear excess of the (3+ + e~ in the energy range of 
~ 300 - 800 GeV. Such a feature was also marginally observed by PPB-BETS^^. 
Furthermore, H.E.S.S. reported a significant steepening of the electron spectrum 
above ^ 1 TcV^"'""''^. The Fermi-LAT Collaboration later reported a high-precision 
measurement of the e+ -|- e~ spectrum from 20 GeV to 1 TeV^^. As shown in 
Fig.l, their measured e+ -|- e~ spectrum may be fitted by a single power-law 
J = (175.40 ± 6.09)e-(3-045±o.oo8) QeV"! m'^ s"! sr"!. Some wiggles are evident 
if one ignores the systematic uncertainties. The spectrum reveals a hardening at 
around 100 GeV and a steepening above ~ 400 GeV. More specifically, the spec- 
trum can be fitted by a broken power-law with indices —3.070 ± 0.025 for £ < 100 
GeV, -2.986±0.031 for 100 < 5 < 400 GcV, and -3.266±0.116 for 400 < 5 < 1000 
GeV^*. So in general, all detectors reveal an excess beyond 100 GeV on the other- 
wise softening spectrum, although different instrumental teams do not agree on the 
detailed spectral shape. 

The other independent indication of the presence of a possible deviation from 
the standard picture came from the recent measurements of the positron-to-electron 
ratio, between 1.5 and 100 GeV by the PAMELA satellite experiment. PAMELA 
found that the positron fraction changes slope at around 10 GeV and begins to 
increase steadily up to 100 GeV^^'^^, confirming the tendency revealed in some ear- 
lier experiments with lower significance and in a narrower energy range (see Section 
2.1). This rising behavior is very different from that predicted for the secondary 
positrons produced by collisions of CR nuclides and the ISM (see cq.{7))^^'^^ . 

We summarize the data in Fig.l and Fig. 2. The latest observations call for an 
additional component of electrons and positrons which is clearly unaccounted for 
in the standard CR model. The lack of an anti-proton excess found by PAMELA^^ 
plays a key role in understanding the origin of the electron/positron excess. 

3. Physics of CRs: Propagation Models and Cosmic Sources 
3.1. The background and the propagation effect 

Galactic cosmic-ray propagation is a decades-old problem in astrophysics. The 
charged particles propagate diffusively in the Galaxy due to scattering with the 
random magnetic field'^'^®'^*. Interactions with the interstellar medium (ISM) and 
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Fig. 1. Observational data against the background model estimates for the e"* 
spectrum (from Ref.26). J is the sum of the fluxes of electrons and positrons. 
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Fig. 2. Some measurements of the positron fraction up to an energy ~ 100 GeV. The data are 
from Refs. 32, 49, 20, 36, 37, 43, 22. The solid line is the positron fraction estimated by using 
eq.(7). 
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the interstellar radiation field would lead to energy losses of CRs. For heavy nuclei 
and unstable nuclei, there are also fragmentation processes by collisions with ISM 
and radiactivc decays, respectively. The distribution of CRs is also modified by the 
convection driven by the galactic wind, and by the re-acceleration due to the inter- 
stellar shocks. As a result, the spectrum of the particles detected on Earth differs 
significantly from those emitted from the sources. The CR propagation equation for 
a particular particle species can be written in the general form^^'^* 

_|[^^_|(v.y)*] (1) 

where ^{f,p, t) is the CR density per unit total particle momentum p at the position 

r, which can be expressed as ^{r,p,t)dp = Airp"^ f{p)dp in terms of phase-space 
density f{p); q{f,p,t) is the source term including primary, spallation, and decay 
contributions; D^x is the spatial diffusion coefficient; V is the convection velocity; 
Dpp is the momentum diffusion coefficient that describes diffusive reacceleration in 
terms of diffusion in the momentum space; p is the momentum gain or loss rate; 
Tf is the timescale for loss by fragmentation; and is the timescale for radioactive 
decay. 

For primary particles such as protons and some heavy nuclei, the source function 
consists of two parts, the spatial distribution /(r) and the energy spectrum q{p). The 
former (/(r)) may follow the distribution of the possible sources of CRs, for example, 
the supernova remnants (SNRs). The injection spectrum q{p) is usually assumed 
to be a power law (or a broken power law) function with respect to momentum p. 
The source function of the secondary particles, which depends on the distributions 
of primary CRs and the properties of the ISM, is given by ^^'^^ 

q{f,p,t) = l3c'^p{f,p,t)[(Tii{p)nuir) + au^nn^if)], (2) 

where /3c is the velocity of the injecting CRs (c is the speed of light); 'i/p{r,p, t) is the 
density per unit momentum p of the primary CRs; cth and cthb are the cross sections 
of producing the secondary particles from the interactions between primaries and 
the H and He targets; nn and riHe are the interstellar hydrogen and helium number 
densities, respectively. The secondary electrons/positrons are mainly produced in 
the proton — proton and proton — He collisions, resulting in charged pions and 
kaons, which further decay as — > tt^ + tt", — > /i^ + u^, — > ji^ + and 
fi^ -\-Vfj^ + Vg. The spatial diffusion coefficient is usually regarded as isotropic 

and is described by a function 

Dxx = pDo{p/po)\ (3) 

where p is the magnetic rigidity (defined as momentum p per unit charge). 

In addition to spatial diffusion, scattering of CR particles on randomly mov- 
ing MHD waves leads to stochastic acceleration, which is described in a transport 
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equation as diffusion in momentum space with a diffusion coefficient Dpp, which is 
related with the spatial diffusion coefficient Dxx through^^ 

D D = ^ (4) 

Here va is the Alfven' speed, and w (may be taken as 1) is the ratio of magnetohy- 
drodynamic wave energy density to the magnetic field energy density, characterizing 
the level of turbulence^°'^^. 

The convection velocity V, which is related to the galactic wind, is assumed 
to be cylindrically symmetric and increases linearly with height z from the galac- 
tic plane"''''. Besides transporting particles, convection also causes adiabatic energy 
losses of CRs as the wind speed increases away from the disk^^. 

Practically most our knowledge of CR propagation is learned from studying sec- 
ondary CRs. This is because the secondary production functions can be computed 
with reasonable precision from the locally observed primary spectra, cross sections, 
and interstellar gas densities. The propagation parameters are determined by com- 
paring the predicted secondaries after the propagation effect with the observations. 
The typical parameters found in the CR modeling are Dq ^ (3 — 5) x lO^^cm^ 
and (5 ~ 1/3 for po ~ 1 GV^^. 

In some simplified cases the propagation equation (1) can be solved analyt- 
ically using the Green function method®^'^^. In most realistic cases an analytical 
solution is not available. A numerical model, widely known as the GALactic PROP- 
agation (GALPROP; http://galprop.stanford.edu) model, was developed by Strong 
and Moskalenko to solve the problem ^9,14 -pj^g main parameters for a given GAL- 
PROP model include the CR primary injection spectra, the spatial distribution of 
CR sources, the size of the propagation region, the spatial and momentum diffusion 
coefficients and their dependencies on particle rigidity. The ISM and interstellar 
radiation field are adopted to calculate fragmentations and energy losses of CRs. 
The parameters are tuned to reproduce the CR spectra observed on Earth. The 
publicly available GALPROP code can give relatively good descriptions of almost 
all kinds of CRs, including the secondaries such as anti-proton as well as diffuse 7 
rays ' . 

As an example, for the Moskalenko & Strong's model 08-005 without 
reacceleration^*, the calculated primary electron and secondary electron and 
positron background fluxes can be parameterized as foUows^^ 

^''S'P™ = ^ TTT(GeV-i cm-2 s'^ sr"!) , 

$bg,sec CGeV"^ cm-2 s"^ sr"^) 

A c:=-0.7 

In the above the energy e is in units of GeV (i.e., e = £/l GeV). For e ^ 1, we 
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have 

$ = $bg,prim ^bg.sec ^bg.sec _ ^bg.prim 3.25 (q\ 

^ e~ e~ e+ e~ ' ^ ' 

and 

-Jlhs = ^bg,sec^^^bg,prim ^ $bg,sec-, _ 0.06£-°-25. (7) 

This set of background spectra agree weh with the results of some more sophisticated 
numerical models^^ derived for e > 10. 

The decreasing behavior of TZ^^ can be understood more straightforwardly. In 
the diffusion models, the diffusion coefficient can be approximated as D oc e^. The 
CRs are assumed to be produced with a power-law spectrum, dN/de oc e~°''=^ . The 
observed spectrum is then a convolution of the source spectrum and propagation 
losses, giving $^5'^"^™ oc £-("e+'5) for primary electrons, where ae and ap are the 
source power-law indices of CR electrons and protons, respectively. Positrons are 
secondary CRs formed from CR protons, and suffer additional propagation loses, 
implying $|^s,scc ^ g-(ap+2<5)_ ^jj^j^ approximately one has TZ^^ cx i<;"<=^"p^'5 as long 
as ^l^^''"^™ > $1^1'^°^ + <j)^s,f=cc^ Usually one adopts a^, so the standard model 
then predicts a CR positron-to-electron ratio that decreases with energy^^'^^'^^, in 
agreement with the tendency reported in eq.(7). 

In the GALPROP code, there are several main approximations: a) The distri- 
bution of sources is continuous, b) All CRs are assumed to be from one single type 
of astrophysical sources, namely the supernova remnants (SNRs). c) The propaga- 
tion of CRs is isotropic. Considering anisotropic propagation models can help to 
solve the problem of soft gradient in the radial dependence of the 7— ray flux, and 
can explain the large bulge-over-disk ratio in positron annihilation as observed by 
INTEGRAL. However, as shown in Fig.4 of Ref.71, the anisotropic propagation 
model gives a similar CR electron spectrum as the isotropic model in the energy 
range > 10 GeV. We therefore will not discuss the anisotropic propagation effect 
in the rest of the review. The change of the first two assumptions can give rise to 
interesting CR spectral signals, which we will come back to discuss later. 

3.2. Sources of cosmic rays 

In general, there are five types of high energy electron/positron production pro- 
cesses: 

• Acceleration of electrons in shocks or in magnetic turbulence. In strong 
shocks the highest electron energy is governed by the balance between en- 
ergy gain (via shock acceleration) and loss (via synchrotron radiation in 
the magnetic field) and can be estimated as fe,max ~ 20 TeV /JghT-Bsh"^^^, 
where Bgh is the strength of the magnetic field in the shock region, /3sh 
is the velocity of the shock in units of speed of light c and F is the bulk 
Lorentz factor of the shock region. 
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• Cosmic rays accelerated at shocks produce secondary electrons/positrons 
inside the source through hadronic interactions. These interactions usu- 
ally lead to production of charged pions or kaons, which decay to pro- 
duce electrons and positrons. Besides the proton — proton and proton — He 
interactions discussed following eq.(3), some other processes include in- 
teractions between hadrons and photons, e.g. proton -|- Tbg — )• A+ — )• 
neutron + 7r+ — > neutron + e+ + v^. + v^, + (proton — 7 process) and 
particle + 7bg — >■ particle + e+ + e^^, for which the energy thresholds are 
fprotoni^Tbg ~ (0.3 GcV)^ and £'particiof7bg ~ WemparticieC^ , respectively. 

• Annihilation of high energy photons with background photons, i.e., photon- 
photon pair production. Thigh + 7bg ^ e+ -|- e~. The energy threshold is 



• Pair production in ultra-strong magnetic fields, i.e., photon-magnetic field 

pair production, 7 + i? ^ e+ + c^. The magnetic field B should satisfy B > 
6 X 10^'^ G (f-y/1 McV)^^ sin??^^. where d is the angle between the photon 
and the magnetic field line^°. Such a strong magnetic field is expected in 
the vicinity of pulsars. 

• Annihilation or decay of dark matter particles. Different from the above 
four processes, such kind of scenarios, although well motivated, are highly 
speculative and lack of solid experimental evidence. 



3.3. Age and distance of the sources of excess 

The above five processes can take place in various astrophysical sources. However, 
it is difficult to pin down the exact source(s) of the observed excess based on 
observations. This is because before reaching Earth the charged particles have been 
deflected by the interstellar magnetic fields. One can see this effect by estimating 
the Larmor radius i?L ~ 10^^ cm {£/! TeV)(i?iG/3/uG)~^, which is much smaller 
than the distance of the CR sources to Earth, where Big ~ 3/uG is the strength 
of the interstellar magnetic field, and Ucmb- Nonetheless, one can reliably estimate 
the age and the distance of these sources in the following way. For an electron with 
energy E^., the total power of synchrotron and inverse Compton radiation in the 
Thomson regime is = io{£ / me(?)^ , where ^0 = 4crTCUtot/3. The characteristic 
cooling timescale can be estimated as ^ f /Pr, i-e.. 
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,2\2 69 
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''rad 



4crTWtot^ 
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where Utot = mb + itcmb + "dust + "star, ub = -Big/Stt is the magnetic field energy 
density, Udust and the photon energy densities of cosmological microwave 

background, dust emission, the star emission, respectively. The travel distance of 
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these CR electrons/positrons cannot exceed the characteristic CR diffusion radius 



~ * ^ 1028 s-1 ^ ' 4 cV cm-3 ^ ' 4o GcV^ ' " ^ ' 

For (5 ~ 1/3, we find that for the observed electrons/positrons at ~ 1 TeV, the source 
should have an age < Trad ^ 10^^ s and should he within a distance < Ras ~ 1 kpc. 
As a result, statistical fluctuations in the injection spectrum and spatial distribu- 
tion of nearby sources within Raa may cause significant deviations from the naive 
predictions in the conventional homogeneous and steady state scenario^^'^^'^^. Since 
i?diff is anti-corrclated to £, the deviations are more prominent at high energies. 

In the following three sections, we will discuss three types of models that have 
been proposed to interpret the observed excess, namely, the modified SNR models 
for the background, new astrophysical sources, and dark matter. 



4. Scenario I: Modified SNR Models for the Background 

Supernova remnants (SNRs) are the canonical sources of CRs'^s. The kinetic energy 
of the SNRs is found to be comparable to that needed to accelerate Galactic CRs. 
The more direct evidence is the detection of the photons up to > 10 TeV from some 
SNRs, which are consistent with 7r° decay from hadronic interactions^^. SNRs have 
been traditionally regarded as the dominant source to contribute to the observed 
background. This model can naturally account for a smooth, gradually softening 
spectrum at high energies. The excess feature discussed above is not expected 
from the simplest model. In this section we discuss whether the excess is simply due 
to some factors that have been ignored in the conventional SNR models. 



4.1. Inhomogeneity of SNRs? 

As already mentioned, in the widely adopted code GALPROP the distribution of 
SNRs is assumed to be homogenous. Such an approximation would not affect the 
results of ion CRs significantly. For electrons, on the other hand, due to the signifi- 
cant cooling caused by magnetic fields and background photons, only the electrons 
generated within a distance ^ 1 kpc can reach Earth'^^'^^ (see also Section 3.3). 
Within such a small volume, the SNR source distribution can be very inhomoge- 
neous. Releasing the homogeneity assumption could therefore lead to significant 
modification of the detected spectrum. 

In the Milky Way, star formation is concentrated in spiral arms^^. The nearest 
spiral arm, the Sagittarious-Carina arm, is about 1 kpc away. The inhomogeneity 
effect of SR source distribution in such a distance scale is therefore very important 
to shape the observed electron/positron CR spectrum^^. Since i?difi is energy de- 
pendent (eq.(9)), an energy-dependent feature would appear. Low energy electrons 

1 GeV) can easily reach us from a larger radius where inhomogeneity is more 
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E [GeV] 

Fig. 3. Predictions of the inhomogeneous SNR models as compared with data (from Ref.66). 
Top Panel: The expected electron and positron spectra - Primary arm electrons (long dashed 
purple), primary disk electrons with nearby sources excluded (short dashed green), nearby SNRs 
(dot-dashed black), secondary positrons (dot-dashed red), and their sum (blue). The hatched 
region describes the solar modulation range (0.2-1.2 GeV). Bottom Panel: Model results and the 
measured PAMELA data for the positron fraction. The shaded region is variability expected from 
solar modulation effects (from Ref.66). 

averaged out. On the other hand, electrons with higher energies (> 10 GeV) can 
arrive Earth only from very nearby sources since they would otherwise have cooled 
via synchrotron and inverse- Compton radiation before reaching Earth. These very 
nearby sources would present a bump in the spectrum in the > 100 GeV range, 
as observed by ATIC and Fermi. Pairs formed in the local vicinity through the 
proton/ISM interactions can reach the solar system also at high energies. As shown 
in Fig. 3, both the PAMELA-observed increase of the positron-to-elcctron ratio be- 
tween 1 and 100 GeV and the ATIC-observed excess around ^ 600 GcV can be 
attributed to the inhomogeneity of SNRs 66,67 _ 

One interesting prediction of the inhomogeneous SNR model is that the positron- 
to-electron ratio should start to drop at energy 8 ^ 100 GeV, just above the present 
PAMELA measurement (see Fig. 3). It should reach a minimum around the "ATIC 
peak" , where it should start to rise again. This is because the nearby SNRs mostly 
contribute to primary electrons near the "ATIC peak" , without significantly increas- 
ing secondary positrons. Such a prediction is in stark contrast to the case where the 
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> 100 GeV excess is due to a primary source of pairs (e.g. either the pulsar model 
or the dark matter models), in which the positron fraction is expected to keep rising 
at a few hundreds GeV until reaching ~ 50%. The positron fraction behavior above 
100 GeV is the smoking gun to test this inhomogeneous SNR model (and the pair 
production models). 

4.2. Accelerution of the secondary in SNRs? 

In the standard model, primary cosmic rays are accelerated in SNRs while the 
secondary are not. In reality the secondary production also takes place in the 

same region where cosmic rays are being accelerated. These secondary participate 
in the acceleration process and turn out to have a flatter spectrum relative to 
the primaries at high energies. This could be responsible, after propagation in the 
Galaxy, for the observed PAMELA anomaly As found in solving the transport 
equation, the downstream equilibrium spectrum of secondary electrons produced in 
the acceleration region has two terms. One term traces the spectrum of the high 
energy primary CRs, i.e., oc £~'^. The other term has a spectrum oc £~'^~^^, where 
the factor arises from the momentum dependence of the diffusion coefficient™'*'^. 
As long as the second term dominates, the spectrum gets flattened. This effect 
was noted earlier as a general expectation for the secondary-to-primary ratio in 
the presence of continuous Fermi acceleration*^'®^. By simulating the spatial and 
temporal distributions of SNRs in the Galaxy according to their known distribution 
statistics, the e+ + e~ energy spectrum measured by Fermi LAT can be interpreted 
within such a secondary acceleration model^° (see Fig. 4). A large portion of the 
model parameter space seems to lead to ovcrprediction of TcV flux as compared 
with the H.E.S.S Observations. Considering the large systematic uncertainties of 
H.E.S.S. and large uncertainties of the model parameters, the model can make a 
case to account for both the excess and the positron fraction excess. Future high 
quality TcV data can give a more reliable test to this model. 

The secondary particle acceleration model has several interesting predictions. In 
particular, all other secondary particles (including baryonic particles) are also accel- 
erated. As a result, the antiproton to proton ratio should also increase with energy. 
This is a clean criterion to differentiate this model from some other astrophysical 
models for positron fraction excess (e.g. pulsar models). Some dark matter models 
however share the same prediction, and therefore cannot be differentiated based 
on this criterion. A similar effect applies to secondary nuclei such as titanium and 
boron These predictions can be soon tested with the data from PAMELA and the 
forthcoming AMS-2 mission. Finally, since this model invokes strong hadronic inter- 
actions, one would naturally expect strong hadronic neutrino signals from nearby 
SNRs, so that IceCube has good prospects to detect them. 

This model requires that secondary positrons can escape the shock region to 
reach Earth. For young SNRs, magnetic field amplification by the shock wave is 
very effective, so that synchrotron cooling of is significant. High energy may 
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Fig. 4. The spectrum of the total cosmic ray electrons and positrons expected in the secondary 
particle acceleration model (The figure is from Ref.80). 

not be able to escape from the acceleration site. As a result, the SNRs invoked in 
this model should be old enough. 

4.3. pair production in the photon-proton interaction in young 
SNRs? 

Despite the cooling constraint discussed above, young SNRs have been argued to be 
an attractive candidate to solve the excess puzzle. In young SNRs, interactions of 
newly formed CR nuclei and background radiation may produce energetic pairs, 
which would change the CR spectrum at the source. Considering a photon field 
at the source peaking in the optical/infrared band {S-y-^,^ ~ eV), pair production 
occurs at (»7icC^/^^7bg)''T'particicC^, which is ~ 1 PeV for protons and ~ a few 
PeV for helium^^. These correspond to the "knee" of the cosmic CR spectrum. In 
the rest frame of the CR nuclei, the secondary electrons/positrons have energies 
around MeV, which turn out to be about TeV in the observer's frame. It is striking 
to note that this is just the energy range where the excesses were observed by 
ATIC, Fermi, and H.E.S.S. Furthermore, the energy density of the excess electrons 
observed by ATIC is about Ugi ~ 3 x 10~^ eV cm~'^ between 0.1 and 1 TeV, which 
is of the order of the energy density due to CR energy loss assuming a spectral 
break from 2.7 to 3.1 at energy ~ 1 PeV. These coincidences are the motivations of 
the photon-CR interaction model*^. Indeed the interaction between CR nuclei and 
background photons may naturally bridge the "knee" of the CR spectrum and the 
excess, provided that the optical background photons have a thermal distribution 
with a temperature Topt ~ 5000 — 7000 K^^ . Such a dense background photon field 
is possible in the vicinity of young SNRs. 
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One potential problem not addressed by Ref.86 is whether the resulting pairs 
can escape from the source. As discussed above, young SNRs tend to have strong 
magnetic fields. Together with the very dense background photons invoked to pro- 
duce e^, strong cooling of the produced secondary is expected, which would 
lead to the production of strong infrared/optical synchrotron radiation and GeV- 
TeV inverse Compton radiation. One needs to check whether the escaped can 
still reach the desired energy to account for the observed excess. 



4.4. Klein-Nishina suppression of electron/positron cooling? 

An electron or positron moving in a dense soft photon bath would lose energy via 
inverse Compton (IC) scattering^^ and produce a high-energy IC component at 

^ic ~ — T+sf*"' "where g = '■ye^^^/me(? and 7o = Elm^-c? is the Lorentz factor of 
the electron. In the Thompson regime, .g ^ 1, so E\c ~ In the extreme 

Klein-Nishina regime, g ^ 1, one has Sic ~ '^e'lTieC^- The scattering cross section 
in the Klein-Nishina regime is also greatly reduced, so that the electron/positron 
cooling rate drops correspondingly. The strength of the drop is proportional to the 
energy density of the seed photons. One may define the critical Klein-Nishina energy 
of electrons £kn ~ 0.3mlc'^/£hg, above which is the Klein-Nishina effect becomes 
significant^*. This corresponds to ~ 44 GeV for a typical B star and ~ 180 GeV 
for typical G-K stars. So given the galactic diffuse optical background from stars, 
the effect takes place in the energy range where electron/positron excess is observed 
to appear^*'^^. Associated with the drop in energy loss rate in the Klein-Nishina 
regime is a gradual hardening of the cosmic ray spectrum. The spectral structure 
in the FERMI and H.E.S.S. data may be reproduced in this way**'*^ provided that: 
(1) The energy density of the seed photons is significantly higher than that of the 
magnetic field (so that IC cooling is more significant than synchrotron cooling) ; (2) 
the star-to-dust emission ratio Wstar/udust is large enough*'^, so that the dominant 
IC cooling is with respect to star optical emission. This requirement is needed since 
IC off the dust emission in the far-infrared range is still in the Thomson regime 
even for relatively energetic electrons. Such a cooling effect would compete with the 
Klein-Nishina cooling effect if Wdust is not suppressed, which would smear the 
spectral feature. This Klein-Nishina-related feature was not noticed previously in 
the GALPROP calculations^^ since a small Ustar/wdust was adopted. 

The Klein-Nishina effect affects both primary electrons and secondary positrons. 
Besides accounting for the excess in the total spectrum, it can affect the positron- 
to-electron ratio and may be able to reproduce the PAMELA result^^. The Wgtar and 
the number density of the interstellar medium are however required to be about 100 
times larger than the average values expected in the Galactic disk. Such extreme 
parameters may apply to the vicinities of young supernova remnants. A caveat is 
again high-energy may not be able to overcome strong cooling and escape from 
such sources (see §4.3). 
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5. Scenario II: New Astrophysical Sources 

In the previous section, it is assumed that SNRs are the sole source of cosmic ray 
electrons/positrons. Such an assumption may well be wrong since other astrophysi- 
cal objects can also produce high energy electron/positron cosmic rays. As early as 
in 1970s, pulsars had been suggested as the source of the CR electrons/positrons^^. 
Some widely discussed objects include Vela and Geminga^"'^^'^^'^^. Some other 
nearby cosmic ray accelerators, such as microquasars or an ancient gamma-ray 
burst, are also possible candidates. In this section we discuss these alternative as- 
trophysical sources. 

5.1. Pulsars 

Pulsars are strongly magnetized, rapidly rotating neutron stars that are powered 
by spin down and emit broad-band piilscid electromagnetic radiation^^. They are 
important particle accelerators. Direct evidence comes from the recent Fermi dis- 
covery of nearly 50 high confidence gamma-ray pulsars with emission above 100 
MeV^*'^^. Three accelerator sites have been discussed in the literature. 1. Polar cap 
models invoke a charge-depleted gap near the magnetic polar cap region^^'^^'^^ in 
which primary particles are accelerated and high energy 7-rays are radiated via 
curvature radiation of inverse Compton scattering. The high energy 7-rays inter- 
act with strong magnetic fields near the pole and produce electron-positron pairs. 
These secondary pairs cool via synchrotron radiation or resonant inverse Compton 
scattering and produce higher generation 7-rays and pairs, leading to a photon-pair 
cascade^^'^'^'^. 2. Outer gap models^'^^'^'^^'^'^^ invoke a charge depleted gap developed 
beyond the "null" charge surface but enclosed by the last open field line. Particles 
are accelerated to high energies in the intense electric field parallel to the magnetic 
field lines inside the outer gap, and 7-rays are produced via inverse Compton scat- 
tering or curvature radiation. The 7-rays interact with other softer photons in the 
magnetosphere to produce pairs. 3. Slot gap models^°^'^'^^ invoke an elongated 
gap near the last open field line, the electric field inside which is immune from pair 
screening thanks to its favorable geometry. Particles can be accelerated from the 
surface all the way to the light cylinder. Almost the entire electric potential can be 
utiHzed for acceleration. Outer gaps cannot exist in old pulsars (e.g. t > 10'' yr), 
while polar gaps and slot gaps can exist in all active pulsars. Recent Fermi observa- 
tions disfavor the polar cap origin of most 7-rays, although both the outer gap and 
the slot gap models are still reasonable candidates for pulsar 7-ray emission, and 
hence, for energetic particle acceleration. Outside the magnetosphere, pairs can be 
re-accelerated in shocks as the pulsar wind interacts with the ambient medium. For 
young pulsar systems, a bright pulsar wind nebula would exist inside the supernova 
remnant, which radiate in a broad band from X-rays up to TeV energies. 

To estimate the CR electron/positron production from pulsars, one needs to 
know the total energy budget of a pulsar and the fraction of this energy that is 
converted into pairs. For the first purpose we adopt the classical rotating magnetic 
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dipole model io6,io7,io8^ -^i^ich the pulsar loses its rotational energy through 
magnetic dipole radiation. The spin-down power is estimated as 

Ldip = mCt ~ 2.6 X 10^° erg s'^ Bl^^Rlf^nl^il + ^/ro)-^ (10) 

where Bp is the polar cap dipolar magnetic field strength of the pulsar, Rg is 
the pulsar radius, fio is the initial angular frequency of rotation (the current 
angular frequency of the rotation is = 27r/P = f2o(l + t/ro)^^/^), tq = 
1.6 X 10^°B~^2^o,3-^45-^s^6 initial spin-down timescale of the pulsar, 

I ^ 10^^ g cm^ is the typical moment of inertia of the pulsar'^'^^. Here the con- 
vention Qx = has been adopted in cgs units. The polar cap dipolar magnetic 
field strength at the surface is Bp = 6.4 x lO^^y/ PP Gauss. Following Refs.109, 
110 we denote /e± as the fraction of the rotational energy that is deposited in the 
cosmic ray pairs. The energy output rate of pairs can be then estimated as 
i?e±,out = fe±I^^- It is usually suggested that the mature pulsars (with ages ~ 10^ 
years and without associated SNRs) rather than the young pulsars can contribute 
to the observed pair flux because the produccid pairs are no longer trapped and 
lose energy in the nebula. The total energy output of the pairs for a mature 
pulsar can be estimated as 

where xq ~ 10^ years/rg. The parameter /g± may be time-dependent and is not well 
determined. A quantitative discussion of plausible values was recently presented in 
Rcf. 109 (Some more model-dependent estimates of the pair production rate 
can be found in Refs.lll, 112, 113). In the context of a standard model of pulsar 
wind nebulae, a reasonable range for /e± falls between 1% and 30%. For typical 
parameters (145, fio, a;o) ^ (1, 300, 10), we have £^e,out ^ 5 x 10^'^(/e± /O.l) ergs. 
To account for the observed ATIC electron/positron excess by a single pulsar, the 
pulsar needs to be within a radius r ~ 0.4 kpc (i^e, out/5 x 10'*'' erg)^/^(ue±/3 x 
10~^ eV cm~^)~^/^. Considering that only a fraction of electrons/positrons can 
reach an energy above 100 GcV, the source should be even closer. Indeed there are 
a few nearby intermediate age pulsars within 0.3 kpc (e.g. Geminga at 0.16 kpc, 
PSR 0656-1-14 at 0.29 kpc, and Vela at 0.29 kpc). For a complete pulsar catalog 
with distance information, see http://www.atnf.csiro.au/research/pulsar/psrcat/. 

Once the total energy output is specified, the injected electron-positron spec- 
trum is fully defined by its spectral shape. The pairs generated from various gaps 
usually are not energetic enough, and they usually do not form a simple power 
law spectrum. In order to be relevant to the CR excess, one needs to assume 
that these pairs are re-accelerated in pulsar wind nebula shocks. It is well known 
that the escaping pulsar wind cannot be accelerated to arbitrarily high energies. 
Depending on the pulsar environment, a cut-off is expected, usually in the energy 
range around TeV. The precise position of the cutoff 7e±,max is rather uncertain, 
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and it critically depends on the pulsar age and on the magnetic field strength in the 
pulsar wind nebula^ 

The pulsar pair spectrum is usually assumed to take the form 

OC 7-±''e-Te± /7e± ,„a. . (12) 

"7e± 

For a point source with a power-law spectrum oc 7^", the energy distribution 
of the pairs reaching Earth can be conveniently handled as^^^ 

ciiVe±/rf7e± [r/Jt;diff(^,7e±)]^e-['-/^^'"(*-^e±)]^ 

«(^'*'Te±) = ^^372;:3 (l-4i7e±)^-« ' ^ ^ 

where the spectrum has a cut off at 7c = I /{tot). The energy loss rate and age then 
set a maximal energy of particles that reach Earth today, with a diffusion radius^ 

RMt,l.±) « 2{D(7e±)ci[l - (1 - 7e±/7c)'-V[(l - ^)7e±/7c]}'/', (14) 

where I'(7e±) « r'o(l + 7e±/7*)'^ and 7, « 6 x lO^. 

Following the above procedure, several authors have calculated the contribution 
of one or more nearby pulsars to the observed spectrum^°^'^^°'^^^. By picking 
a few nearby known pulsars and employing some suggested energy output models, 
one is able to reproduce the spectral features and the intensities of the reported 
excess, with reasonable values for the output efficiency (see Fig. 5 for illustration). 
The advantage of this model is the lack of an anti-proton excess, in agreement with 
the PAMELA observation^''. Since electrons and positrons are produced in pairs 
in the pulsar models, the positron fraction would achieve 50% at high energies 
in the pulsar model. Another prediction of the pulsar model for the excess 
is that at higher energies, the electron/positron flux measured at Earth will be 
dominated by a few nearby pulsars, and therefore the spectrum would host wiggle- 
like features^^'^°^'^^°. The presence of such features at high energies would strongly 
suggest a pulsar origin of the anomalous contribution to the electron/positron fluxes. 
The hardening of the; CR electron spectrum detected by H.E.S.S at an energy ~ 4 
TeV, if confirmed in the future, may be taken as a support to the pulsar model, 
although the inhomogeneous distribution of SNRs may also produce the wiggle-like 
structure in the electron/positron spectrum at TeV energies. The differentiation 
could be made by evaluating the corresponding positron-to-electron ratio, which 
should be small for the SNR model (since secondary positrons are much rarer than 
primary electrons), but large for the pulsar models (electrons and positrons are 
produced in pairs). 

5.2. Microqusars 

Microquasars are a sub-group of accretion-powered X-ray binaries that process a 
relativistic jet. The companion can be either a high mass or a low mass star. The 
accretor itself is likely a rapidly spinning black hole, although a neutron star accretor 
is not ruled out. These objects display in miniature some of the properties of quasars. 
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Fig. 5. The spectrum of cosmic ray electrons and positrons expected in the pulsar model (from 
Ref.109). 



accreting supermassive black holes with much higher luminosities, and hence, carry 
the name of "microquasar"^^''. 

Although their physical origin is still unclear^^®, microquasars can be important 
sources of high energy cosmic rays^^^'^^'^. Gamma-ray flares with energy above 100 
MeV have been detected from microquasars LS I +61 303^^^ and Cygnus X-3^^^, 
suggesting that they are powerful lepton accelerators. TeV electrons/positrons may 
be accelerated in the internal or termination shocks. More pairs can be produced as 
TeV gamma-rays attenuate with the UV/optical photons from the companion star. 
The net energy output of microquasars in the Galaxy is uncertain. An optimistic 
estimate suggests that up to 5 — 10% of the total CR luminosity can be contributed 
by such kind of objects^^^'^^°. If correct, the energy output of microquasars is high 
enough to account for the detected CR excess, as noted in Refs.24, 21. However, 
there are not many known nearby microquasars^^'^. The closest candidate is Cen 
X-4 at ^ 1.2 kpc. There is no detailed calculation of this model in the literature. 



5.3. Gamma ray Bursts 

Gamma-ray bursts (GRBs) are brief, intense flashes of soft (0.01 — 1 MeV) 7- 
rays in the sky, which are related to collapses of some massive stars or mergers 
of two compact stellar objects (two neutron stars or a neutron star and a black 
hole) They were serendipitously discovered by Vela-satellites in late 1960s^^"', 
and now are routinely detected by gamma-ray detectors and followed up by vari- 
ous telescopes in broad band (from radio to TeV). The event rate of nearby high- 
luminosity GRBs (those typical GRBs detected at cosmological distances) is roughly 
0.5 — IGpc^'^yr"^, which corresponds to about 0.025-0.05 GRB per Myr per 
galaxy^^^. Counting for those GRBs that do not beam towards us (with large un- 
certainty in the beaming correction factor) - which are relevant since off-beam 
GRBs can also contribute to the diffuse background - the rate can be as high as 
1 GRB per 10^ — 10^ yr per galaxy. Nearby low-luminosity GRBs can have a higher 
observed event rate (by a factor of 200-500)^^^'^^°. Even if their collimation angles 



August 30, 2010 0:10 WSPC/INSTRUCTION FILE electron-positron 



22 Y. Z. Fan , B. Zhang & J. Chang 

are typically wider, the coUimation-corrected event rate could be still higher than 
that of high- luminosity GRBs. As a result, it is possible that an ancient GRB in the 
solar neighborhood has produced enough in a not distant past, which diffuse to 
Earth, leading to the observed cxcess^^^. Such a model can interpret both the total 

excess and the positron fraction excess. Furthermore, depending model param- 
eters, both a sharp feature measured by ATIC and a smooth feature measured by 
Fermi can be rcproduced'^'^^ . 

There are however two caveats related to this scenario. First, it is not clear 
whether our galaxy is suitable to host cosmological GRBs. Observations show 
that cosmological GRBs typically reside in host galaxies that are different from 
p^j.gi32,i33 gome augthors argued that Milky Way is too metal rich for bright 
GRBs^^**, while some recent observations indicate that at least some GRBs can re- 
side in host galaxies with metallicity as high as 0.3 solar value^^^. More data are 
needed to settle down the debate. Second, even if bright GRBs can occur in our 
galaxy, it is non-trivial to deposit nearly same amount of 7-ray energy to TeV 
pairs that escape from the GRB. One may envisage several pair production mecha- 
nisms in a GRB: (i) A baryonic fireball may be pair-rich in the early stage. However, 
most pairs would have annihilated when the fireball reaches the 7-ray radiation ra- 
dius, (ii) In the prompt 7-ray emission region, TeV photons would annihilate with 
hard X-ray photons to produce TeV pairs, but these secondary pairs would quickly 
cool and re-radiate in the MeV regime; (iii) Prompt MeV photons may annihilate 
with the soft gamma-rays back-scattered by the circum-burst medium to produce 
pairs. These pairs have much lower energies. They may be reaccelerated by the ex- 
ternal shock, but the escaping fraction may not be high, (iv) TeV photons emitted 
by the GRB (likely from external shock since the compactness parameter is too 
high for internal prompt emission) may interact with the "background" UV/optical 
photons to produce TeV pairs. These pairs can most easily escape. The question is 
how to get "background" UV/optical photons. Prompt bright optical flash similar 
to that detected in the naked-eye burst GRB 080319B^'^^ was invoked in the model 
of Ref.131. Such bursts are very rare, so that the chance of having a nearby GRB 
capable of accounting for the observed excess is further reduced. 

6. Scenario III: New Physics 

Although the electron/positron excess can be explained within the astrophysical 

scenarios, another widely discussed scenario is that the excess is the signature of 
dark matter. The dark matter models have strong astrophysically motivations, and 
can interpret the excess data under the assumption of some (extreme) parameter 
regimes. In this section, we critically review these models. We begin with a bricif re- 
view of the observational evidence of dark matter, and some widely-discussed dark 
matter candidates (see also Refs. 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 
148, 149 for extensive reviews on dark matter candidates). We then go through 
various dark matter models (both annihilation and decay models) proposed in the 
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literature to interpret the electron/positron excesses, and comment on the observa- 
tional constraints to the parameter spaces of these models. For a detailed review on 
the annihilation models, see also Ref.27. 

6.1. Observational evidence for dark matter 

It was recognized as early as 1930s^^° that some matter in the Universe is invisible. 
The non-baryonic nature of dark matter was established since early 1970s^^^. In 
the following we summarize several compelling arguments for the existence of non- 
baryonic dark matter. 

6.1.1. Galactic rotation curves 

The most direct evidence for dark matter in the galactic scale comes from the obser- 
vations of the rotation curves of spiral galaxies (i.e., the graph of circular velocities 
Vc of stars and s function of their distance r from the galactic center), which 

are usually obtained by measuring Doppler shifts of the 21 cm emission line of hydro- 
gen gas. The detected rotation curves usually exhibit a characteristic flat behavior 
at a distance that is beyond the edge of the visible disks (see Fig. 6 for illustration). 
This implies that the mass M<(r) enclosed within a radius r follows Af<(r) oc r 
for r ^ rdisk, suggesting the existence of an extended dark matter halo. Measuring 
stellar velocity dispersion in elliptical galaxies also leads to a similar conclusion. 
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Fig. 6. Measured rotation curve of NGC6503 with best fit and the contributions from the dark 
matter halo, disk and gas (from Ref.153). 



August 30, 2010 0:10 WSPC/INSTRUCTION FILE electron-positron 



24 Y. Z. Fan , B. Zhang & J. Chang 

6.1.2. Galaxy Clusters 

Galaxy clusters are the largest gravitationally bound objects in the Universe. Three 
independent methods all suggest that dark matter is the main mass that bind 
clusters together. 

First, Newton's gravity law suggests that the velocity dispersion of galaxies is 
approximately v^{r) ~ GM^(r)/r. One can infer the total mass M<(r) from the 
measurements of the velocity dispersion and the size of a cluster (which can be 
determined if the rcdshift and angular size of the cluster arc measured). In 1933, 
Zwicky^^° used this method to measure the mass-to-light ratio of the Coma cluster 
to be 400 times of the solar value. This gave the first hints of dark matter in the 
modern sense. The conclusion is confirmed by many follow-up measurements. 

Second, according to the virial theorem, the gas temperature of the clusters 
may be estimated as k^T « 1.5 keV (M<(r)/10^'' Mq){1 Mpc/r), where M<(r) 
is normalized to the baryonic mass derived from the star light M{, 10^^ Mq. 
However, the detected temperature is as high as ~ 10 keV, suggesting the existence 
of much more massive dark matter, i.e., M<(r) ~ 6Mb. 

Finally, mass can be measured via the gravitational lensing effect predicted 
by Einstein's general theory of relativity^^*. The mass of a lensing cluster can be 
inferred by the deflection angle 6 « (GAfci/rfc^)^''^ and the impact parameter d, 
both can be measured or inferred if the redshift of the lensing cluster is known. 
Usually the derived M^x is much larger than the observed baryonic mass M^, again 
indicating the c;xistence of the large amount of invisible matter. 

The most strong support to the dark matter pictiire came recently from an 
X-ray observation of the bullet cluster (Fig. 7). The X-ray image (red) shows the 
distribution of baryonic mass (in the form of hot gas). The bulk of mass distri- 
bution as inferred by the gravitational lensing technique is denoted in blue, which 
does not trace the distribution of hot gas. Such a configuration is difficult to inter- 
pret within other competing models (e.g. the modified Newtonian dynamics), and 
strongly supports the dark matter interpretation. 

6.1.3. Cosmological constraints on energy densities 

Advances in technology have usher in the new age of precision cosmology. Wilkinson 
Microwave Anisotropy Probe (WMAP) can precisely measure the primordial fiuc- 
tuation in the temperature map of the cosmic microwave background (CMB). By 
precisely locating the first acoustic peak in the power density fluctuation, WMAP 
clearly measured the curvature of the unvierse^^^. Together with the observational 
results of Type la supernova that suggest an accelerating universe^^'', very precise 
values of the parameters of the flat ACDM have been obtained. According to the lat- 
est 7- year WMAP data^^*, the energy densities of various components of the ACDM 
model are fl^ = 0.0456 ± 0.0016, = 0.227 ± 0.014, and f^A = 0.728 ± 0.015. 
An independent constraint on = fidm + f^b from the baryon acoustic peak is 
consistent with the WMAP result^^^. All these strongly suggest that the universe 
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Fig. 7. The X-ray image by Chandra X-ray Observatory of the bullet cluster (red) over-imposed 
on the mass distribution (blue) derived from gravitational lensing. It is clearly shown that 
most of the mass does not trace the gas emission in the X-ray band. This Figure is from 
http: / /www. spacoimages.com/hubuclle06ph. html. 

is dark-matter dominated, and that the bulk of dark matter is not of the baryonic 
origin. An independent constraint on the amount of barynic contribution to dark 
matter was derived from the MAssive Compact Halo Objects (MACHO) experi- 
ment, which showed that the micro- lensing detected MACHOs can only account for 
a small fraction 20%) of the missing mass in the Milky Way and LMC halos^^*'. 

6.2. Dark matter properties and possible candidates 

The evidence for dark matter is compelling in all astrophysical scales. To be a 
viable dark matter candidate, one needs to satisfy some general requirements^"*^: 
(1) Dark matter must be dark, in the sense that it must have no (or extremely 
weak) interactions with photons. This is effectively saying that the electric charge 
of dark matter particles must be extremely smaller than the electron charge. (2) 
Self-interactions of the dark matter should be rather weak. This is manifested in the 
dark matter map of the Bullet Cluster (Fig. 7). The collision between two galaxy 
clusters cause strong shocks of baryonic gas between the two clusters, while the 
two dark matter halos streamed through each other without noticeable interaction. 
(3) Interactions between dark matter and baryons must be very weak. Otherwise 
baryons and dark matter would fall together in the overdense region and form a 
baryon-DM disk in galaxies rather than forming an extended DM halo. (4) Dark 
matter cannot be mainly made up of standard model particles, since most leptons 
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and baryons are charged. The only potential candidate is neutrinos, but they are 
too "hot" and too light to be trapped in DM halos. 

Physicists arc never short of ideas of inventing non-baryonic dark matter 
candidatcs^^^'^""^. Almost all current models of dark matter use the standard concept 
of quantum field theory to describe the properties of elementary particle candidates 
(for exceptions, see for instance Refs.161, 162). This means that the dark matter 
particles can be characterized by their mass and spin. In the following we summa- 
rize the properties of several widely discussed dark matter candidates. As shown in 
Table 1, the mass of the proposed candidates spans a wide range. 

The most studied dark matter candidates are weakly-interacting massive par- 
ticles (WIMPs). The main motivation for introducing WIMPs is the gauge hi- 
erarchy problem, namely, the observational constraint of the mass of the spin 
Higgs Boson (a hypothetical massive particle in the Standard Model (SM) in par- 
ticle physics), i.e. 114.4 GeV < WHiggs < 186GeV^^'''-^^'*, is much lighter than the 
Planck mass Mpi = 1.2 x lO^^GeV. This suggests that SM is not the complete 
model, and there should be new physics in the energy range of the weak inter- 
action scale mweak ~ 10 GeV — TeV. WIMPs are the predicted particles of the 
new physics, which naturally carries a mass with the similar energy scale. Very 
interestingly, this energy scale is broadly consistent with the energy scale of the 
PAMELA/ATIC/HESS/Fermi excess. This was the main reason that the 
excess stimulated the great excitement from the particle physics community. Three 
types of WIMPs, i.e. ncutralinos and other supersymmetric particles, Kaluza-Klein 
particles in universal extra dimensions, and inert Higgs doublets, are discussed be- 
low. Some candidates, such as axions and sterile neutrinos, have too low an energy 
to be relevant to the excess data. Since they arc physically well motivated, we 
nonetheless include them for the sake of completeness. 

Tabic 1 . Properties of some Dark Matter Candidates (adapted 
from Refs.147, 149). 



Type Particle Spin Approximate Mass Scale 



Neutralino 


1/2 


10 GeV - 


10 


TcV 


Kaluza-Klein UED 


1 


TeV 






Inert Higgs Doublet 





~50 GeV 


or 


~650 GeV 


Axion 





;tieV-meV 






Sterile Neutrino 


1/2 


keV 







6.2.1. Neutralinos and other supersymmetric particles 

The most widely discussed WIMPs are neutralinos, one type of the supersymmetric 
(SYSY) particles. Supersymmetry (SUSY) is an ingredient in many superstring the- 
ories trying to unite all the fundamental forces of nature, including gravity The 
gauge hierarchy problem is most elegantly solved by SUSY. In SUSY extensions of 
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the standard model, every standard model particle has a new, as yet undiscovered 
partner particle, which has the same quantum numbers and gauge interactions, but 
differs in spin by 1/2. From the standpoint of dark matter, the lightest supersym- 
metric particle is naturally stable in models that conserve i?-parity, which is defined 
as 

i?= (-1)3^5+^+2^, (15) 

where B is the baryon number, L the Icpton number and S the spin of the 
particle^'^*'^^^'^'*^. The ordinary particles and their superpartners have R = +1 
and —1 respectively, which means that supcrsymmetric particles can only be cre- 
ated or annihilated in pairs in reactions of ordinary particles, or a single supcr- 
symmetric particle can only decay into final states containing an odd number of 
supersymmetric particles (plus some Standard Model particles). The idea that su- 
pcrsymmetric particles could be good dark matter candidates became attractive 
when it was realized that breaking of SUSY could be related to the electroweak 
scale, and that the supersymmetric partner of the photon (the photino) would cou- 
ple to fermions with electroweak strength^^^. Below we list some supersymmetric 
dark matter candidates^^'''^^^ with brief discussion. 

• Neutralinos in models of R-parity conserving supersymmetry are by far the 

most widely studied dark matter candidates. Since the superpartners of 
the Z boson (zino), the photon (photino) and the neutral higgs (higgsino) 
have the same quantum numbers, they can mix to form four eigenstates of 
the mass operator called "neutralinos" In many models the lightest of 
the four neutralinos x turns out to be the lightest supcrsymmetric particle 
(LSP), which is, therefore, stable and can only be destroyed via pair an- 
nihilation (Neutralinos have spin 1/2 and are their own antiparticles; that 
is, they are Majorana fermions). x can be produced thermally in the hot 
early imiversc and leave approximately the right relic abundance to account 
for the observed f2dm/i^ ~ 0.1, where h ~ 0.71 is the Hubble's constant in 
units of 100 km s~^ Mpc~^. These facts make x with a mass of roughly 
10-10000 GeV, an excellent cold dark matter candidate. They may be de- 
tected directly through scattering in detectors, or indirectly through the 
decay products that result when neutralinos annihilate in pairs^^^'^^^. 

• Sneutrinos are the superpartners of the Standard Model neutrinos in the 
siipersymmetric models. They have a spin 0, and have long been considered 
as dark matter candidates'^''. If their mass is in the range of 0.55 — 2.3 TeV, 
the sneutrinos can have a cosmologically interesting relic density'^®. From 
the direct detection searches sneutrinos have been excluded as the major 
component of the dark matter'^^. However, in the supersymmetric models 
with the inclusion of right-handed fields, sneutrinos are still compatible 
with the current direct detection sensitivities and could be a viable dark 
matter candidate'^°. 
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• Gravitinos are the spin 3/2 supcrpartners of gravitons in the supersym- 
metric models, and could be the lightest, stable supersymmetric particle. 
Although strongly theoretically motivated, gravitinos as dark matter would 
be very difficult to observe due to their extremely weak interactions with 
ordinary matter^^^. Long lived gravitinos may impose problems to cosmol- 
ogy. For example, their presence may destroy the abundances of primor- 
dial light elements^''^'^'''^. These problems can be circumvented in some 
scenarios^''^'^^'^. For instance, in the case of a broken R-parity, the con- 
straints from primordial nucleosynthesis are naturally satisfied and decay- 
ing gravitinos would lead to characteristic signatures in high energy cosmic 

174 

rays ' . 

6.2.2. Kaluza- Klein Particles in Universal Extra Dimensions 

Our world appears to consist of 3 + 1 dimensions (three spatial dimension and one 
time dimension). It is however possible that there exist other 6 spatial dimensions, 
which appear at much higher energy scales, as suggested by Kaluza and Klein in 
the 1920s. In the extra-dimension models, the (3 + l)-dimcnsion space-time one 
experiences is a structure called a "brane" , which is embedded in a (3 -|- J -h 1) 
space-time called the "bulk" . Upon compactification of the extra dimensions, all the 
fields propagating in the bulk have their momenta quantized in units of ^ 1/=K^, 
where K denotes the size of the extra dimcnsions^^^. As a result, for each bulk field 
a set of Kaluza-Klein (KK) states (i.e. Fourier expanded modes) appear. In the 
3-1-1 dimensional world, these KK states appear as a series of states with masses 
TTin = where n labels the mode number. Each of these new states contains 

the same quantum numbers, such as charge, color, etc^'*^. Usually the standard 
model fields are assumed to be confined on the brane and only gravity is allowed 
to propagate in the bulk. However it might be possible for all the standard model 
particles freely propagate in the extra dimensions^°°. Such a scenario is the so-called 
universal extra dimensions. The simplest universal extra dimension models preserve 
a discrete parity known as the KK-parity, implying that the lightest KK particle 
(LKP) is stable and can be an interesting, possible dark matter candidate^'*^'^''^. 
The LKP is typically the level 1 partner of the hypercharge gauge bosons, denoted 
as 200 Tjjg WMAP limit of fldm/i^ = 0.1131 ±0.00341^8 corresponds to a mass 
of the dark matter candidate B^^^ between roughly 0.5 to 1 TeV, depending on the 
exact form of the mass spectrum and the resulting co-annihilation channels^^^^^^^. 
Current collider measurements give a constraint of > 0.3 TeV, whereas LHC 
may probe the KK mass parameter space up to 1.5 TeV^^^. 

6.2.3. Inert Higgs 

Inert Higgs are a type of WIMPs that are introduced in one of the most minimal 
extensions of the standard model. Back in 1978, it was noted that a model with two 
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Higgs doublets containing a discrete symmetry could contain a state, the lightest 
neutral scalar or pseuodoscalar boson, which is stable^''^'^^^. It is the subject of the 
renewed interest recently since it could be used to solve the naturalness problem 
in the standard model^^^ . The possibility of one of the lighter neutral states in the 
enlarged Higgs sector to be the dark matter was also pointed out, and soon the basic 
properties of this "inert" Higgs candidate for dark matter were investigated^^°. It 
turns out that this model contains a dark matter candidate, an "inert" particle 
that couples gauge bosons only and is stable due to its discrete symmetry. Taking 
into account the relic density constraints from WMAP and various theoretical and 
experimental constraints, the latest constraint on the mass of the Inert Higgs is 
~ 30 — 80 GeV (middle mass range) or ^ 500 — 800 GeV (high mass range) 
Candidates in the middle mass range may annihilate through Z or through Higgs. 
Some solutions in the middle mass range may also have a large annihilation cross 
section, through loop corrections, into a pair of gamma rays. This is relevant for 
indirect detections of dark matter through annihilation at the Galactic Center, and 
such a gamma ray line would be easily observed by the Fermi satellite^^^. For the 
inert Higgs dark matter candidates in the high mass range, annihilation into 
pairs is allowed with a large cross section^^°. 

6.2.4. Axions 

Axions were introduced in an attempt to solve the strong CP problem in particle 
physics^^^. They are still one of the leading dark-matter candidates. The axion mass 
is defined by 



where TOu — 4 MeV, md ~ 8 McV and w 135 McV arc masses of the up quark, 
down quark and pion, and /jr « 93 MeV is the pion decay constant. The axion 
phenomenology is therefore determined by, up to a numerical factor, by one num- 
ber only, i.e. the energy scale /a of symmetry breaking. A variety of astrophysical 
observations and laboratory experiments constrain /a to be between 10^ GeV and 
10^^ GeV, suggesting an axion mass Wa ~ 10~^ — 10~^ eV. Smaller masses would 
lead to an unacceptably large cosmological abundance, while larger masses are ruled 
out by combinating of constraints from supernova 1987A, globular clusters, labora- 
tory experiments, and the search for two-photon decays of relic axions^^^'^''^. The 
lifetime of axion is T(axion -> 77) « ^^^^-p — where is a model-dependent 

parameter and 5^ is expected not to deviate from 1 significantly. In general, its very 
small mass makes axion not possible to interpret the 100 GeV excess discussed 
in this review. 

Within the SUSY theories, axions have spin 1/2 superpartners, which are called 
axinos. Models combining SUSY and the Peccei-Quinn^''^ solution to the strong 
CP problem necessarily contain this particle. Depending on the model and the 
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(16) 
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SUSY breaking scheme, the axino mass can range from eV to GeV. Although widely 
believed to be only capable of acting as a warm, or hot dark matter candidate^^^'^^*, 
it was recently shown that for quite low reheating temperatures, cold axino dark 
matter may be possible^*'^'^*^. Interesting astrophysical signals are expected from 
axinos if the R-parity is broken^^''. 



6.2.5. Sterile Neutrinos 

The Standard (active) neutrinos were considered as a favored particle dark matter 
candidate in late 1970's, when the calculations of the massive neutrino relic density 
were first carried out^^*'^^^'^^°'^^^. To account for the dark matter of a dwarf galaxy 
of velocity dispersion C7d typically of order 100 km s~^ and core radius Tc typically 
1 kpc, the neutrinos have to be "massive" i93,i46 

.n. > 120 eV f l^^f ri^y (17) 



which is much larger than the upper limit m,^ < 2.3 eV derived from the various 
experimental results (e.g. the tritium decay experiments at Mainz^'^^. Therefore 
the standard neutrinos are too light to be the dominant component of dark matter, 
although among the proposed candidates for non-baryonic dark matter standard 
neutrinos have the "undisputed virtue of being known to exist" 

Sterile neutrinos, unlike standard (active) neutrinos, do not interact through 
standard weak interactions^^'''^^^. Instead, they communicate with the rest of the 
neutrino sector through fermion mixing^^®'^^''. The existence of sterile neutrinos 
(right-handed or gauge singlet) is one of the most attractive explanations of the ob- 
served flavor oscillations of active neutrinos. Sterile neutrinos were proposed as dark 
matter candidates in 1994^^^^. Stringent cosmologieal and astrophysical constraints 
on sterile neutrinos come from the analysis of their cosmologieal abundance and 
the study of their decay products^^^. A very recent analysis finds out that sterile 
neutrino DM with mass > 2 keV is consistent with all existing constraints"'^®^. On 
the other hand, this mass is too low to account for the 100 GeV excess feature 
discussed in this review. 



6.3. Dark matter models for the excess 

In many dark matter models, the dark matter particles are their own anti-particles 
(Majorana fermions), and can annihilate with each other. There are several possi- 
ble DM annihilation scenarios: (i) DM particles annihilate to standard model par- 
ticle pairs, such as gauge bosons, quarks, and lepton pairs^^^'^^^'^°^'^°^'^°^'^''^; (ii) 
DM particles annihilate via virtual internal bremsstrahlung processes to produce 
e"'"e~7^°^; (iii) DM particles annihilate to new mediating particle pairs, which then 
decay to e+e~ and other standard model particles^'^^'^°^'^°®. According to these 
scenarios, dark matter annihilation would lead to electron/positron pairs, and may 
account for the PAMELA/ATIC/HESS/Fermi e^ excess (see e.g., Refs.210, 211, 
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212, 138, 213, 141, 27, 207, 214, 215, 216, 217, 218, 206, 219, 220, 221, 222, 223, 
224, 225, 226, 227, 228, 229, 230, 231, 232, 233, 234, 235, 236, 237, 238, 239, 240, 
241, 242, 243, 204. 244, 245, 246, 247, 248, 209, 249, 250, 251, 252, 253, 254, 255, 
256, 257, 258, 259, 260, 261, 262, 263, 264, 265, 266, 267, 268, 269, 270, 271, 272, 
273, 274, 275, 276, 277, 278, 279). According to these scenarios, the resulting lepton 
spectra would have a cut off energy at the DM mass mdm- Hence, the cutoff energy 
seen in Fcrmi/HESS (or ATIC) data would suggest mdm 1-5 (or 0.6) TcV. 

The dark matter self-annihilation scenarios, which have been extensively dis- 
cussed since 1970s, are not the only possibilities to produce signatures for dark 
matter indirect detections. Strictly speaking, the viability of a particle as a dark 
matter candidate does not require its absolute stability, but merely require that 
the dark matter lifetime Tdec is much longer than the age of the Universe (i.e., 
Tdec ^ 10^^ s). Hence, if dark matter decays at a sufficiently high rate, the decay 
products might be detectable. In particular, if the DM decays predominately into 
leptons in a time scale Tdec ~ 10^^ s (which is long enough not to cause other 
cosmological problems), the detected 100 GeV excess may be also reproduced 
(see e.g., Refs.173, 280, 61, 281, 282, 283, 284, 285, 286, 287, 288, 289, 290, 291, 
292, 293, 294, 295, 296, 297, 298, 299, 300, 301, 302, 303, 304, 305, 306, 307, 308, 
309, 310, 311, 312, 313, 314, 315). In these DM decay scenarios, the resulting 
spectra would have a cutoff at an energy ^ mdm/2. Hence, the cutoff energy seen 
in Fermi/HESS (or ATIC) data would suggest mdm 3 (or 1.2) TcV. 

In both the annihilation and decay scenarios, the DM rcquircid in the modeling 
belongs to the WIMP category. This is encouraging since weak-scale particles make 
an excellent dark matter candidate, having the merit of the so-called WIMP miracle: 
In the simplest type of models of thermally produced dark matter, an average 
of the annihilation rate at the time of chemical decoupling can be estimated as 
< av 3 X 10~^^(0.1/f2dm^^) cm^ s~^, where < av > is the thermal average of 
the total annihilation cross section a multiplied by velocity v. Taking the relative 
velocity of colliding DM particles at the decoupling time v ^ 10^° cm s~^, one has 
a ~ 10-8 GeV-2 ~ "V^dm' mdm ~ 100[a/(l/137)] GeV. This mass is in the 
energy range of weak interaction energy scale, and hence the particle is likely a 
WIMP. 

In the following, we review the annihilation and decay models in turn. 



6.3.1. Dark matter annihilation models 

In the dark matter annihilation scenarios, the source term of the electrons/positrons 
is given by 

,^ X l/pdm(r)\^ dN^± 

where is the energy spectrum of the electrons/positrons produced by DM an- 
nihilation, pdm(r) is the dark matter density distribution in the annihilating region, 
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and < av >~ 3 x 10~^^ cm'^ s~^ based on the WIMPs miracle argument discussed 
above (i.e. the usual DM annihilation rate producing the excess signal is also 
related to the annihilation rate producing the cosmological relic DM density). With 
a given q{f,p), one can solve tho propagation equation (i.e., Eq.[l]) and get the 
final electron/positron spectra detected on Earth. For such a purpose, the density 
distribution of the dark matter halo is needed. Various models of DM halo density 
distribution have been proposed. Some popular ones can be parameterized as 



where pdm,© ~ 0.3 GeV cm is the DM density in the vicinity of the sun, ^ 8.5 
kpc is the solar distance from the Galactic center, and the profile parameters 
(w, C, Tf", rs) are (2, 2, 0, 5 kpc), (1, 3, 1, 20 kpc), (1, 3, 1.5, 30 kpc) for 
the Core Isothermal model^"*^^, the Navarro, Frenk and White (NEW) modeP^^ and 
the Moore model^"'^^, respectively. It is well known that only the TeV CR pro- 
duced within a distance ~ 1 kpc can reach us, i.e., r << r©, hence the above 3 
popular (but different) dark matter density distribution models give rise to rather 
similar CR signals. The photon signals from these three different models, such as 
the GeV- TeV gamma-rays generated from inverse Compton scattering off the CMB 
photons by the TeV e^'s or the synchrotron radiation of the TeV e^'s, can how- 
ever reach us from a much larger distance range including the core of the DM halo. 
These signals therefore can be very different for different models. If DM annihila- 
tion also produces hadronic pairs, the spectrum of the antiproton CRs can be also 
very different for different models, since antiprotons are more rigid and much less 
radiative so that they can reach from much larger distances. The non-detections 
of a prominent electromagnetic spectral signals in the radio (synchrotron origin) 
and GeV- TeV (inverse Compton origin) bands, as well as the non-detection of a 
prominent antiproton cosmic ray signal, therefore impose stringent constraints on 
the DM density profile and/or the DM annihilation properties. 

Since the first release of the PAMELA data, many candidate WIMP DM models 
have been studied. The spectral shape of the observed excess is not difficult 
to reproduce if the annihilations proceed mostly to leptons rather than to quarks, 
gauge bosons, and so on. Eor example, two-body annihilations directly into pairs 
produce a hard spectrum, which can be adjusted to be consistent with the ATIC 
data. If the two-body final states are /U^ or pairs, their subsequent decays into 

can produce a soft enough spectrum consistent with the Eermi and HESS data. 

However, in all these models, the thermal average cross section multiplied by 
velocity < av >= 3 x 10~^^ cm^ s~^ is found to be too small to account for the flux 
level detected by PAMELA/ATIC/Ecrmi/HESS. An unnaturally large "enhance- 
ment/boost factor", Ep ~ 10^ — 10^, is needed. Quite a few mechanisms have been 
proposed to introduce such a Ep factor. (1) If DM halos have substructures, the local 
density can be enhanced so that the annihilation optical depth can be increased^^^. 
Detailed N-body simulations, on the other hand, show that the boost factor intro- 




(19) 
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duced by DM substructures is generally less than Ep ~ 10; (2) It is possible that 
DM in the early universe was not in equilibrium. In the scenarios of non-thermal 
production of DM, the interaction rate responsible for the excesses are not di- 
rectly related to the relic DM density. One may simply assume that the interaction 
rate is that required from the excess data, and no boost factor is needed^^"; (3) 
More physically, there might exist new mediating particles during the DM annihi- 
lations. Examples include the Sommerfeld mechanism^^-'^'^^^ and the Breit-Wigner 
mechanism^^^. The Sommerfeld mechanism^^^ is a non-relativistic quantum me- 
chanical effect. A large boost factor can be produced if DM interacts with a new 
light particle. Without such an interaction, annihilation of DM is a short distance 
effect. If DM interacts with a nearly masslcss (light) particle, then there is an ad- 
ditional long-range force (with Coulomb-like potential V{r) = —h/r), which would 
distort the DM wave-function before annihilation happens. This would introduce a 
modification factor to the annihilation cross section, so that a — cr^S. The factor S 
is the Sommerfeld factor. By solving the non-rclativistic Schrodinger equation, one 
can derive = S. Assuming a dark matter particle x coupling to a mediator 
with coupling strength A, for S-wave annihilation in the nonrelativistic limit, one 
gets Ep = (7r5/v)/(l — e~^^/^) for = and Ep < Trbrndm/mip for > 0, where 
V is the velocity (in units of c) of each DM particle in the center-of-mass frame^'^''. 
Ep ~ 10^ is then achievable if v ~ 10~^67r for = or mdm/'m^ ~ lO'^/bn for 

> 0. The Breit-Wigner mechanism^^^ applies to the case of a heavy mediating 
particle. For most masses of the mediating particle, the cross section is not en- 
hanced. However, if the annihilation of DM is through the S-channel and if the new 
particle has a mass about twice of the DM mass, the cross section and annihilation 
rate can be greatly enhanced near a resonance point. This introduces a large enough 
Ep to interpret the data. 

Another challenge the DM model faces in the excess modeling is the non- 
detection of an excess in the antiproton spectrum^^. Contrary to the data, most 
DM annihilation models invoke a large hardronic annihilation fraction. 

The requirement of having excesses only in but not in anti-proton spectrum, 
along with the large boost factor needed as discussed above, pose great constraints 
and eliminate many candidate DM annihilation models. For example, the most 
popular DM candidate, the lightest supersymmetric particle (LSP) in the minimal 
supersymmetric standard model (MSSM) can give rise to the excesses^'^^. How- 
ever, as a linear combination of photino, zino and higgsino, LSP usually has a large 
hadronic annihilation fraction, which is in conflict with the PAMELA anti-proton 
data. This model also suffers the problem of realizing a large boost factor. For 
the KK dark matter particles, direct annihilation to e+e~ is possible. In fact, the 
lightest KK particles (LKP) was recognized as a potentially important source 
of positrons when the excess was found by HEAT earlier. In the ATIC discov- 
ery paper^^, the KK models was also introduced as an example to account for the 
observed e+-|-e~ spectrum peak around 600 GeV. These models again require a sig- 
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nificant fraction of annihilation in hadronic modes (~ 35%), which is at odds with 
the PAMELA antiproton data (cf. Ref.240). The required large boost factor E-p is 
attributed to halo sub-structures, which was found difficult to achieve. In the inert 
Higgs model, for a candidate with mdm ^ 70 GeV, the DM annihilation contribu- 
tion to and anti-proton fluxes is much smaller than the expected backgrounds. 
Even if a boost factor is invoked to enhance the signals, the particle candidate is un- 
able to explain the observed e+ and anti-proton data. Only a high mass candidate, 
with mdm ~ 10 TeV, would be possible to fit the PAMELA excess provided a large 
enhancement, either in the local DM density or through the Sommerfeld effect^^^. 
Such a heavy particle, however, is in contradiction with the constraint from WMAP 
and various theoretical and experimental eonstraints^^^. 

Recently two types of models have been constructed to account for the detection 
of the excess in the electron/positron channel but not in the antiproton channel (see 
Ref.27 for a recent review): 

(1) The; kinematieally limited light particle decay models (e.g., Refs.207, 242, 
243, 204, 244, 245, 246, 247, 248, 209, 249, 250, 251, 252, 253, 254, 255, 256, 257, 258, 
259, 260, 261, 262); These models require the existence of a light particle with a mass 
less than the sum of proton and anti-proton masses. Such a low mass makes hadronic 
modes kinematieally inaccessible and forces the particle decay predominantly into 
leptons. To accommodate the very high annihilation rate required to reproduce the 
current excess, the light particle decay model makes use of its light force-carrier 
to generate the Sommerfeld enhancement^"^. For instance, in the MSSM extended 
by a singlet chiral superfield, the singlino-like neutralino dark matter annihilates to 
light singlet-like Higgs bosons, which proceed to decay to either electron-positron 
or muon-antimuon pairs. For the singlino-like neiitralino dark matter with a mass 
~ 2 TeV, PAMELA and Fermi data may be interpreted^'^^'^^^. Phenomenologically, 
the light particle decay can have multiple channels. For example, for the case of 
the mediating light particle being a light scalar i/)^"^'^"^, for annihilation XX ~^ 4'4': 
one would consider four (/) decay channels: 1) — > e+e~, 2) (/> — >■ fi'^ fi^ , 3) a 
mixture of 1:1 (or 2:1) between electrons and muons, and lastly 4) </> — )• 7r+7r~^°^. 
Detailed modeling suggests that the models annihilating predominantly to fj, are 
somewhat favored^^^. For the case of mediating light particle is a scalar s and 
a pseudoscalar a, i.e. xX ~^ ^c-^ one has a also preferably decay to muons, i.e. 
a A*"*"/^", and muons further decay to detected on the Earth. Some spectrum 
modeling examples^^^ can be found in Fig.8, where the model N3 is one of the 
Nomura-Thaler-type models^*^ and is characterized by = 0.8 GeV and nis = 20 
GeV. The model AH4 is one of the Arkani-Hamed et al.— type models^°^, in which 
the scalar dark matter particles annihilate into n'^ n~ ■ 

(2) The leptophilic DM models2i4.2i5, 216,217,218 xhese models can be realized by 
demanding DM interactions being leptophilic, namely, mediating DM particles only 
interact with the standard model leptons. The mediating particle does not have to 
be light. The only constraint on the mass comes from the requirement of sufficient 
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Fig. 8. Example spectra of good fit DM annihilation models found in Ref.324. The predicted 
e"*" + e~ signal and background are compared against the Fermi and HESS data. The HESS data 
and the background model have been rescaled with a factor 0.85. In the inset, the positron fraction 
as measured by PAMELA is plotted against the predicted signal for the same model. 



enhancement of the DM annihilation cross section. In the model with Sommerfeld 
enhancement the force-carrier particle should have a mass rridm/ E-p, while in the 
model with Brei-Wigner enhancement the force-carrier particle should have a mass 

The DM annihilation models, although most intensively discussed in the litera- 
ture, are now facing strong observational and theoretical constraints. If the required 
boost factor Ep > 10'^ is not due to the DM substructures, and if one assumes such 
a boost and adopts a universal Navarro- Frenk- White (NEW) density profile of DM 
halos, a large annihilation 7— ray flux from nearby galaxy clusters is predicted. With 
the parameters that well fit the PAMELA and Eermi data'^^'* and the standard as- 
sumptions for the limiting mass of the substructures within DM halos, it is found 
that the EGRET upper limit in 7-ray flux is violated in Virgo^^^. Similarly, the pre- 
dicted gamma-ray flux violates the HESS upper limit constraints from the Galactic 
center and some dwarf spheroidals, and the predicted radio flux from synchrotron 
emission of the produced from the Galactic Center also violates the current radio 
observation upper limits. The upper limits may be satisfied only if the DM density 
profile is significantly less steep than the benchmark NEW and Einasto profiles^^^. 
The diffuse 7— ray data from the Eermi first year observations have imposed a very 
tight constraint on the DM annihilation models and suggests that the annihilation 
models likely only work in some "fined-tuned" situations'^^^'^^^. Erom theoretical 
point of view, it is shown^^^ that the enhancement required in the Sommerfeld- 
enhanced DM annihilation models implies a thermal relic density that is too small 
to allow the DM candidate suitable for interpreting the feature to account for 
all the DM in the universe. In Ref. 328 the authors also derived an upper bound on 
the possible Sommerfeld enhancements from the observations of elliptical galactic 
DM halos and showed that these bounds generically exclude such an explanation. 
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Counterarguments were also raised recently. It was argued"^^*^ that there are a lot 
of uncertainties in calculating the 7-ray signals from the nearby dwarf spheroidals 
and the Galactic Center, and that the annihilation models do not conflict with the 

Fermi gamma-ray limits. It was also found that the annihilation models are better 
than other models to account for the microwave signals collected by WMAP in the 
Galactic Center^^^. 



6.3.2. Dark matter decay models 

The dark matter particles are often assumed to be perfectly stable as the result 
of a symmetry, e.g. R-parity in the supersymmetric models. However, from the 
gravitational evidence for the existence of dark matter we can only infer that dark 
matter has to be stable on timescales (much) longer than the age of the Universe. 
In other words, DM particles are allowed to decay. If they decay at a sufHciently 
high rate, the decay products might be detectable as cosmic rays, which would re- 
sult in the observed excess. There are in fact some physically well-motivated 
dark matter candidates that have very long lifetimes before decaying. For instance, 
gravitino dark matter, which is unstable due to a small breaking of the R-parity, 
constitutes an interesting scenario that leads to a thermal history of the universe 
consistent with the observed abundances of primordial elements, the observed dark 
matter relic abundance, and the observed baryon asymmetry Some other de- 
caying DM candidates discussed in the literature include hidden gauge bosons^*^, 
hidden gauginos^®^, right-handed sneutrinos in R-parity breaking scenarios^®^, KK 
dark matter in a simple extension of the minimal universal extra dimension mod- 
els by introducing a small curvaturc^'^^ , and baryonic bound states of messenger 
quarks^*'^. The late decay of dark matter particles, for example gravitinos, can pro- 
duce gamma rays^^°'^°^, positrons neutrinos^^^ and/or antiprotons^^^, which 
contribute to the; total fluxes received at the Earth. 

In the dark matter decay scenario (s), the source term of the electrons/positrons 
is given by 

,{r,p) = ^'^^^'-^, (20) 

Tdec JTldm 

where '^'^^ is the energy spectrum of the electrons/positrons produced by the 
decay of each DM particle. Evidently, in the decay model the parameters involved 
regarding the DM properties are the mass m^^, the lifetime Tdec, and the spectrum 
of decay, 

In order to account for the excess in the electron/positron channel but not 
in the antiproton channel, the decay of DM should mainly produce leptons with 
suppressed hadronic branching ratios. As a result, the leptophilic decaying DM 
scenarios have attracted considerable attention. As in the annihilation scenarios, 
leptophilic decaying DM models can be broadly divided into two categories'^*. One 
is that DM first decays into some new light particles, which subsequently decay 
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into muons or electrons, while decays into hadrons are forbidden by kinematics. 
The other is that the DM particle couples mainly to leptons due to symmetry^^^ or 
geometric setup'^"®. 

Although the decay models are not as widely investigated as the annihilating 
models, they have some advantages that are not shared by the annihilation models. 
First, in the decay models there is no need of the "boost factor" given the free- 
dom of choosing the decay lifetime of the DM particles. Second, decay models suffer 
much less constraints from the gamma-ray and radio observations from the Galactic 
Center^^®'^^^'^^^. Similar to the annihilation models, only the TeV pairs gener- 
ated within a distance 1 kpc can reach us. In the popular NFW profile of the 
dark matter halo with p{r) cx r~^, one expects a much higher DM density near the 
Galactic Genter. The gamma-ray signals produced in inverse Gompton scattering 
off the star light and CMB by pairs produced in the Galactic Center can reach the 
Earth. Since the decay flux is proportional to the DM density p rather than propor- 
tional to as in the annihilation models, the decay models predict a much fainter 
signal than the annihilation models. The non-detection of such signals pose strong 
constraints on the annihilation models, but not the decay models (c.f. Ref.336, 330). 

In the phenomenological approach, it is usually assumed that the dark matter 
particle is either a fermion {ipdm) or a scalar {(pdm)- Predictions for the positron 
fraction for various decay channels and different dark matter masses and lifetimes 
can be then computcid^^^-"*'^^. In the case of a fermionic dark matter particle, people 
usually consider the two-body decay channels ipdm Z^v, '0dm ^ W^l^, as well 
as the three-body decay channels ipdm — >■ with i — e, p, t being the charged 

leptons. For a scalar dark matter particle, one usually considers the two-body decay 
channels (pd^ Z°Z°, 0dm -^W+W', and 0dm ^ In light of the PAMELA 

and Fermi data, some specific channels are found better than the others^"^^. One 
example of spectral modeling is shown in Fig. 9. The latest diffuse 7— ray data from 
Fermi first year observations excludes the dark matter decay model unless the dark 
matter decays mostly or even exclusively in ji^ (Fig. 10, see however Refs. 327, 
338 for looser constraints). 

7. Discussions and Conclusions 

After 40 years of efforts, the cosmic ray electron/positron spectrum has been 
measured up to ^ 1 TeV/ ~ 100 GeV with an unprecedented accuracy. The rising of 
the positron-to-electron ratio above ^ 10 GeV, a tendency discovered in the 1970s^^, 
has been confirmed by the PAMELA satellite^^ with high confidence. The electron 
-|- positron spectrum hardens above a few hundred GeV, as discovered by ATIG^^, 
PPB-BETS^s, Fermi-LAT^s and HESS'"'"'^^ The rising positron-to-electron ratio 
and the hardening electron -|- positron spectrum, unexpected in the standard cosmic 
ray model, are consistent with each other and are called "excesses" . The detailed 
electron -|- positron spectral shape recorded by ATIG/PPB-BETS and Fermi/HESS 
are however different. Future more advanced experimental observations will pin 
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Fig. 9. Positron fraction (left panel) and electron+positron flux {right panel) for the fcrmionic 
DM decay (i/'dm ~^ M^M"'') rnodel. The parameters are mjin ~ 3.5 TeV and t^cc ~ 1-1 x 10^® 
sec. The dashed hnc shows the astrophysical background. A scalar dark matter particle decay 
(0dm M"*"/^" ) model gives an equally well fit to the data for rridm ~ 2.5 TeV and Tdec ~ 1.8 X lO'^^ 
sec. The figures are from Ref.337. 

down the shape of the spectrum more accurately. 

A lot of efforts have been made to interpret the observed eletron/positron ex- 
cesses. We summarize them in three main categories: 

• In the standard cosmic ray model, the supernova remnants are the sole 
source of primary cosmic rays up to PeVs. The inhomogeneity of the su- 
pernova remnants^®'®'' and/or the secondary particle acceleration in su- 
pernova remnants™'®'^ may account for the PAMELA/ATIC/Fermi data. 
More speculatively, some physical processes taking place in young super- 
nova remnants, for example, suppression of the high energy cooling in 
the Klein-Nishina regime^^'*® and pair production by interactions be- 
tween high-energy CRs and background photons®^, may also account for 
the observed excesses. 

• It is possible that the observed excesses originate from a new astro- 
physical CR source. The leading candidate is pulsars^"'^^^^^^. Pulsars are 
known to be able to produce electron/positron pairs, which may be acceler- 
ated to the desired high energies outside the magnetospheres (likely in SNR 
shocks). No excess in the spectrum of antiproton cosmic rays is expected, 
in agreement with the PAMELA result. The escape of these pairs from the 
associated SNR, however, is unclear. One popular assumption is that pul- 
sars should be mature, for which the associated SNR is too old to trap the 
pairs. The potential problem is that mature pulsars have lost a significant 
fraction of their initial kinetic energy and may not be powerful enough to re- 
produce the observations. Some other new astrophysical sources discussed 
in the literature include microqusars^^^ and local Gamma-ray Bursts^'^^. 
The results are rather uncertain in these scenarios. 

• The most extensively discussed interpretation of the excesses is the anni- 
hilations or decays of dark matter. The excesses in positrons and antiprotons 
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Fig. 10. The dark matter model parameter space — Tdcc (dark matter mass— lifetime) and 
the constraints from the diffuse galactic gamma ray measurements by the Fermi satellite in its 
first year operation (from Ref.326). The exclusion contours are due to the Fermi observations of 
the '10° — 20° 'strip', the |6| > 60° 'Galactic Poles' region, and the isotropic flux, respectively. 
The regions allowed to fit the PAMELA positron data (green and yellow bands, 95% and 99.999% 
C.L. regions) and the PAMELA positron + Fermi and HESS data (red and orange blobs, 95% and 
99.999% C.L. regions) in terms of decaying Dark Matter are also plotted. 
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in the CR spectra are two "standard" predictions of these models ' . 
The detection of a rising positron-to-electron ratio is a good news for these 
models. The non-detections of an excess in the antiproton spectrum, how- 
ever, calls for modifications of most canonical dark matter models. Phe- 
nomenologically, one usually assumes that^^^'^^'''^^'^'^^ (1) the dark matter 
particles annihilate or decay into new particles, which are so light that 
can only finally decay into and/or /x^; (2) the DM is leptophilic and 
couples only to leptons due to symmetry. In the DM annihilation mod- 
els an unnatural large "boost factor" ^ 10^ is needed to reproduce the 
PAMELA/ATIC/Fermi detected fluxes. Such an unusually large boost 
factor is difficult to interpret simply by invoking substructures in the DM 
halo. Some particle-physics-oriented suggestions, such as the Sommerfeld 
and Briet-Wigner mechanisms, have been invoked to account for the large 
excess. Many specific models have emerged and are able to fit the elec- 
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tron/positron excesses. There are however some additional observational 
constraints. For example, almost all dark matter models predict a cer- 
tain level excess of 7 ray. The latest diffuse 7— ray data from the Fermi 
first year observations have excluded a large portion of the DM parameter 
space^^^'^^'^'^^^. (In contrast, the constraints on the astrophysical models 
are rather weak, since they are "local" phenomena and just occupy a small 
fraction of the volume of the Galaxy). In the case of annihilating DM, those 
models generating predominantly /x^ and assuming a cored Isothermal-like 
DM profile can fit the data with a mass mam ^ 2 TeV, but for other param- 
eters, the model suffers various observational constraints and is less favored 
(c.f. Ref.330). In the case of decaying DM, again except those models gen- 
erating predominantly /x^ (which can fit the data with mdm ~ 3 TeV and 
Tdec ~ 2 X 10^^ sec), most models are not favored by various observational 
constraints. It is worth noting that recently direct searches of dark mat- 
ter have been performed by various experimental groups. For example, the 
CDMS II collaboration has observed two candidate events"^"^^. If interpreted 
as DM recoil on nuclei, the DM particles detected by CDMS II should be 
relatively light, with a mass mdm ^ 0(100) GeV. This is much lower than 
that needed in the excesses modeling. Due to low statistics, it is of course 
too early to say that the observed excesses are not of the DM origin. In 
any case, direct detection experiments and the high energy colliders (such 
as the Large Hadron Collider) would play an essential role to probe the na- 
ture of DM, and to provide a test to the current DM interpretations of the 
excesses. In the near future, neutrino observations by the upcoming high 
energy neutrino detectors (such as IceCube) would also help to constrain 
the DM nature of the excesses^"^^'^^^. 

All the models reviewed here assume that the detected cosmic ray excesses 
are due to one physical process. In reality, this may not be the case. For example, it 
is possibly that the detected excesses are mainly contributed by astrophysical pro- 
cesses, and the DM annihilation/decay signals are outshone. If so, the extraction 
of the DM properties via indirect detection experiments would bo oven more chal- 
lenging. Nevertheless further indirect detection experiments are still necessary since 
they are complementary to the direct detection experiments and the collider exper- 
iments for revealing the properties of DM. For instance, collider experiments may 
identify a long-lived, weakly interacting particle and measure its mass. On the other 
hand, it would be difficult to test its cosmological abundance and stability. Such 
goals, however, are achievable in the direct and/or indirect detection experiments. 

We conclude that both astrophysical processes and DM annihilation or decay 
models can provide a self-consistent explanation to the recently observed excesses 
in the cosmic ray spectrum. At present the dark matter interpretations seem to be 
more exotic than some astrophysical scenarios. This is in particular the case for 
the dark matter annihilation models which are tightly constrained by the 7— ray 
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observational data by Fermi. In fact the latest diffuse gamma-ray data have excluded 
a large region of the DM parameter space that is favored in interpreting the 
excesses. More experimental observations are needed to pin down the natin-o of DM 
and its implications for the excess. Among the promising astrophysical models, 
the inhomogeneous distribution of supernova remnants has a smoking gun prediction 
that the positron-to-electron ratio is expected to decrease above an energy ~ 100 
GcV. The pulsar model, instead, predicts an increasing positron-to-clcctron ratio 
up to the energy ~ 1 TeV. These two predictions will soon be tested by the ongoing 
PAMELA and the upcoming AMS-2 experiments. Both the inhomogeneous SNR 
model and the pulsar model predict a possible wiggle-like structure in the TeV 
energy range, which may be constrained by the observations of H.E.S.S-like ground 
based Chercnkov telescopes. 
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